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We explore how to extract effective dynamics from loop quantum gravity and spinfoams trun- 
cated to a finite fixed graph, with the hope of modeling symmetry-reduced gravitational systems. 
We particularize our study to the 2-vertex graph with N links. We describe the canonical data 
using the recent formulation of the phase space in terms of spinors, and implement a symmetry- 
reduction to the homogeneous and isotropic sector. From the canonical point of view, we construct 
a consistent Hamiltonian for the model and discuss its relation with Friedmann-Robertson- Walker 
cosmologies. Then, we analyze the dynamics from the spinfoam approach. We compute exactly 
the transition amplitude between initial and final coherent spin networks states with support on 
the 2-vertex graph, for the choice of the simplest two-complex (with a single space-time vertex). 
The transition amplitude verifies an exact differential equation that agrees with the Hamiltonian 
constructed previously. Thus, in our simple setting we clarify the link between the canonical and 
the covariant formalisms. 
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Introduction 



Loop quantum gravity [l| and spinfoams |2| form together a proposition for a well-defined framework for quantum 
gravity. While loop gravity is the canonical definition of the theory describing the evolution of quantum states of 
space geometry, the spinfoam approach provides the covariant formulation of the theory, such that it describes the 
quantum structure of space-time. More precisely, in both loop quantum gravity and spinfoams the quantum states 
of geometry are given by spin network states which have support on some graphs. The space of spin networks living 
on all possible graphs (up to diffcomorphisms) provide a basis for the kinematical Hilbcrt space of the theory. In the 
canonical framework, the evolution is implemented by a Hamiltonian operator. There exist both graph-changing and 
non-graph-changing proposals for this operator depending on the precise regularization scheme and implementation 
of the Hamiltonian constraint at the quantum level, although it is usually assumed that it acts on the states by 
changing the underlying graph. The spinfoam approach, on the other hand, defines transition amplitudes between 
spin network states living on arbitrary graphs through the construction of a covariant discretized path integral. This 
program towards quantum gravity faces three main issues: a clear definitive definition of the dynamics, the derivation 
and analysis of the semi-classical regime of the theory where we should recover fluctuations of the gravitational field 
around flat space-time, and a consistent method for extracting quantum gravity corrections and predictions. 

Here, we propose to discuss these topics in the context of loop quantum gravity and spinfoams truncated to a finite 
fixed graph. From the canonical point of view, this requires defining a Hamiltonian on a fixed graph (without assuming 
that it comes from the truncation of a graph-changing or a non- graph-changing Hamiltonian). From the covariant 
point of view, this requires considering transition amplitudes between initial and final spin network states with support 
on the same graph. The theory restricted to spin network states living on this given graph, is thus truncated to a 
finite number of degrees of freedom. Our hope is that restricting the theory to a finite fixed graph would allow to 
formulate physically relevant mini-supcrspacc models for (loop) quantum gravity. Indeed, mini-supcrspacc models 
in general relativity are restrictions to certain families of 4-metrics parameterized by a finite number of parameters 
and satisfying certain symmetries or properties which made them relevant to some particular physical context. We 
expect that the development of such mini-superspace models of loop quantum gravity will lead to realistic models for 
quantum cosmology and thus allow precise cosmological predictions from loop gravity and spinfoam models. 

Our strategy is to choose some simple graph T, to analyze and describe both the classical phase space and the 
space of quantum states living on this graph, to define and study the dynamics both at the classical and quantum 
levels using the loop gravity ansatz for the Hamiltonian or the transition amplitudes of spinfoam models, and finally 
to understand how it can be mapped (or not) on some cosmological models (or other interesting situations). Working 
on such simplified setting with a fixed underlying graph allows to define rigourously the dynamics of the classical 
data and quantum states and more generally to investigate the possible dynamics that one can define. We also 
hope that studying such toy models for (loop) quantum gravity will allow to understand more about the geometrical 
interpretation of the quantum states and the construction of coherent states, about the transition to the semi-classical 
regime, and about the structure of spinfoam amplitudes for the evolution of spin networks. From this point of view, 
the fact that such models could lead to realistic cosmological models or to the dynamics of other symmetry-reduced 
geometries, and that we could possibly extract quantitative quantum gravity effects in this context would be a bonus. 

We discuss the generic procedures and methods behind this fixed finite graph approach, but we focus in practice 
on the case of the 2-vertex graph, which has been shown to be somewhat related to Friedmann-Robertson- Walker 
cosmology in earlier works • The first section describes in detail the classical kinematical structures of loop gravity 
on an arbitrary fixed graph T. We review the recently developed approach of parameterizing the classical phase space 
with spinor variables [7|, l9Ull| and discuss its relation to the other parameterizations in terms of the standard loop 
gravity holonomy-flux variables, in terms of twisted geometries foL fl"3| . and finally in terms of U(.ZV)-covariant variables 
[g, 0, IT3 - HH . Each set of variables allows to insist on certain aspect of the kinematics and clarifies the geometrical 
interpretation of the phase space. 

In the second section, we apply the generic method to the particular case of the graph with 2 vertices and N edges 
linking them, which is the simplest graph on which one can formulate the theory. We describe the phase space on 
this 2-vertex graph and define the symmetry reduction to the homogeneous and isotropic sector, following the U(iV)- 
symmetry proposal of @, • Then in section three, we define and implement classical dynamics on this 2-vertex graph 
consistent with the reduction to the homogeneous and isotropic sector. We show the relation between this truncated 
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loop gravity classical dynamics, the (geometrical part of the) Hamiltonian for Fricdmann-Robcrtson- Walker (FRW) 
cosmology and the effective dynamics derived from loop quantum cosmology. Let us remark that here we focus on 
the vacuum case. In the future we will need to include matter in order to get true models for cosmology. We further 
discuss and explain the limitations of the similarities between our model and cosmology, and the failures of our naive 
2-vertex graph Hamiltonian at large scales when we take into account a non-vanishing curvature or cosmological 
constant. 

In the fourth and final section, we investigate the effective classical dynamics on the 2-vertex graph induced by the 
quantum transition amplitudes of spinfoam models applied to coherent spin network states on the 2-vertex graph. 
This clarifies and extends the previous results obtained by Bianchi, Rovelli and Vidotto in We compute exactly 
the spinfoam transition amplitudes and identify the exact differential equations that they satisfy. Then we show 
how these differential equations lead back to the 2-vertex Hamiltonian defined previously at the classical level. We 
discuss how to generalize these results beyond the 2-vertex graph and how our procedure is related to the study of 
recursion relations and invariance of spinfoam amplitudes jlmlq . Indeed, from the spinfoam point of view, such 
recursion relations are understood to be equivalent to the dynamics of the theory and have been shown to translate to 
differential equations when applied to coherent spin network states. These are exactly the differential equations that 
we recover in our 2-vertex graph setting and that we show to be related to the Hamiltonian constraint of flat FRW 
cosmology. We also discuss how to modify the spinfoam amplitude to take non-vanishingcurvature into account. This 
final step finally shows the coherence of the spinfoam cosmology approach initiated in [5[ with the canonical point of 
view, although much work is needed to go beyond our naive truncation to the 2-vertex graph and its homogeneous 
and isotropic sector. 

I. CLASSICAL PHASE SPACE OF LOOP GRAVITY 

Loop quantum gravity is formulated in terms of spin network states leaving on graphs. A spin network state on a 
graph r is defined as a gauge-invariant function of SU(2) group elements g e leaving on the edges e 6 T. The SU(2) 
group elements are physically the holonomies of the Ashtekar-Barbero connection along the edges of the graph. The 
Hilbert space of these wave- functions <^r(<7e) provides a quantization of the phase space of holonomy-flux variables 
(g e ,X e ), where the holonomies g e act by multiplication and fluxes X e act as derivation operators. Following recent 
developments on the V(N) formalism for intertwiners 0, 0, fi"3T[l5l EH and twisted geometries [H, H^ . it has been 
understood that the holonomy-flux algebra defined in terms of the (g e ,X e ) variables on a given graph T can be 
re- written using spinor variables z v e living around each vertex v on the edges e 0, lij-fllj . Then wave-functions will be 
holomorphic functions of these spinor variables. 

In this section, we will quickly review this construction. We will define the phase space of the spinors z^, provide 
an action principle encoding the canonical Poisson structure and constraints generating the SU(2) gauge-invariance, 
and explain how to recover the standard holonomy and flux observables from these variables. Finally, we will discuss 
how to endow this kinematical structure with dynamics, thus defining effective classical dynamical models of loop 
quantum gravity on fixed graphs. 

A. Spinor Networks and Phase Space on a Fixed Graph 

In loop quantum gravity the spin network states provide an orthonormal basis for the kinematical Hilbert space 
of the theory. A spin network state has support on a given closed graph and consists in the coloring of the graph's 
edges and vertices with appropriate quantum numbers. More precisely, every edge is colored with an irreducible 
representation of SU(2) and every vertex with an intcrtwiner, namely a SU(2)-invariant state living in the tensor 
product of the SU(2) representations meeting at that vertex. This kinematical Hilbert space is usually formulated as 
the quantization of the phase space of holonomy-flux variables, which are discretized observables for the connection 
and triad field of loop quantum gravity. Recently, new descriptions for the kinematical phase space of loop gravity 
have been devised in order to understand better the geometrical interpretation of the spin network states, also with 
the aim of constructing suitable coherent states of discrete geometry. Such descriptions, as the twisted geometries and 
the U(iV) formalism for intertwiners, have converged to a description of the phase space in terms of spinor networks. 

Let us start with a closed oriented graph T, with E oriented edges and V vertices. For simplicity's sake, we choose it 
connected, else all the definitions will still apply to each connected component of the graph. Now around each vertex 
v, we associate a spinor variable z v e £ C 2 to each edge e attached to v. Equivalently, this amounts to associating to 

each edge e two spinors, zj = z| and z\ = Ze , the former attached to the source vertex s(e) of the edge and the 
latter attached to its target vertex. 
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The phase space is defined by the canonical bracket on the spinor variables, postulating that each spinor z is 
canonically conjugated to its complex conjugate: 

{z a , Zb} = -iSab , (1) 

where we have dropped the e, v indices and z a with a = 0, 1 stand for the two components of the spinor z. Then we 
impose two sets of constraints on these sets of spinors: 

• Closure constraints at each vertex v : 



W, C„=^|<)«|--«K)I = 0, (2) 

where \z)(z\ and I are 2x2 matrices and the linear combination \z)(z\ — is the traceless part of 

• Matching constraints along each edge e: 

Ve, M e = {z s e \z s e )-{z t e \z t e )=0. (3) 

It is direct to check that these two sets of constraints are first class. The closure constraints generate SU(2) 
transformations at each vertex: 

K) — ► g v K) with g v e SU(2) , (4) 

while the matching constraints generate U(l) transformations on each edge: 

\zl) -> e+^|<), |4) e-^\4), with e+*- € U(l) . (5) 

Moreover, one easily checks that closure and matching constraints commute with each other. A spinor network is a 
set of spinors {z^} satisfying both sets of constraints and up to SU(2) and U(l) transformations, i.e an element in 
C 4B //(SU(2) 1/ x XJ(l) E ), where // stands for the symplcctic reduction (i.e both solving the constraints and quotienting 
by their action). 

This constrained phase space structure can be summarized by an action principle: 

si%z] = f dt^^(-z(z:i^:) + (z:iA„|z:)) + ^p e ((z:iz:)-(z*i^)), (6) 

v e3v e 

where the 2x2 matrices A.„ with Tr A„ = and the real variables p e are Lagrange multipliers imposing the closure 
and matching constraints. This action defines the kinematical structure of spinor networks and the phase space 
on the fixed graph T. We will later add a Hamiltonian term to this action in order to define (classical) dynamics for 
these spinor networks, the goal being to construct the Hamiltonian in order to produce effective dynamics for loop 
quantum gravity on fixed graphs T which can be relevant for symmetry-reduced physical situations such as cosmology. 

In order to understand the geometrical meaning of spinor networks, the best is to translate spinors into 3- vectors. 
Indeed, each spinor z € C 2 determines a 3-vector FgI 3 through its projection onto the Pauli matrices: 

V=(z\a\z), \z)(z\ = \((z\z)l + V-a), (7) 

where the Pauli matrices are normalized such that of = I, Vi and the norm of the 3-vector is |V| = (z\z). Reversely, 
the spinor z is entirely determined by its projection V up to a global phase (see in appendix for more details). 
Swapping all the spinors z v e for their projections V£ , it is straightforward to translate the constraints in terms of the 
3- vectors: 

• Closure constraints at each vertex v. 



(8) 
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• Matching constraints along each edge e: 

Ve, \V?\ = \V:\ (9) 

The geometrical interpretation then appears clearly. The closure constraints mean that each vertex v is dual to a 
(closed) polyhedron with each edge e dual to a face of the polyhedron. The 3-vector V" becomes the normal vector to 
the face dual to e and the norm | V" | gives the area of that face. Every set of vector satisfying the closure constraint 
automatically defines a unique such polyhedron. A detailed reconstruction of the polyhedron from the normal vectors 
is achieved through Minkowski theorem and Lasserre's algorithm (2pj . 

Thus we have one polyhedron (embedded in the flat 3d Euclidean space M 3 ) around each vertex v. Then the 
matching constraints impose that the areas of their matching faces along each edge e be equal. Note that this does 



not mean that the shape of the faces will match. This would require further constraints (see e.g. |21|). This translation 
of the closure and matching constraints in terms of 3-vectors provides spinor networks with a clear interpretation as 
discrete 3d geometries. 

Now, we would like to define obscrvables and wave-functions over the classical spinor phase space C 4B / /(SU(2) x 
U(1) B ). Since we have two sets of constraints, there are two natural paths to solving them depending if we first 
implement the U(l)-invariance or the SU(2)-invariance. To start with, on the initial unconstrained phase space C 4B , 
we have wave-functions ip(z^). which can be defined simply as holomorphic functions of the spinor variables. Then 
we have two alternatives: 

• We first impose the matching constraints A4 e : 

This is the path to the standard formulation of loop (quantum) gravity and to twisted geometries. Natural 
U(l)-invariant observables, constructed from the spinors, are group elements g e £ SU(2) attached to each edge 
e. Together with the 3-vectors V£, they allow to entirely parameterize the phase space C 4E / /XJ(l) E . These g e 
define the SU(2) holonomies along edges of loop gravity and can be taken as the configuration space coordinates. 
We will review the reconstruction of these group elements from the spinors in the next subsection If B[ 

Then we would work with wave-functions f(g e ), which are U(l)-invariant. Finally imposing the closure con- 
straints on these wave-functions, we would obtain a Hilbcrt space of SU(2)-invariant wave-functions f(g e )-, with 
the standard spin network functionals as a basis. 

This scheme seems to localize the degrees of freedom (at the kincmatical level) on the edges of the graph. 

• We first impose the closure constraints C v : 

This is the path taken by the U(7V) formalism for SU(2) intertwiners @, Ell, EH ■ Imposing the closure constraints 
at each vertex v, one defines natural SU(2)-invariant observables F^, depending on pair of edges e, / attached 
to v, and holomorphic in the spinor variables z" e . Their definition will be reviewed below in subsection IIP I 

Then we would work with wave- functions ip(F^), which are SU(2)-invariant. Finally imposing the matching 
constraints on these wave-functions, we obtain a Hilbcrt space of U(l)-invariant wave-functions <p(F£f) on the 
graph r, which is exactly isomorphic to the standard Hilbert space of spin networks obtained by the other path 
of first imposing U(l)-invariance and then SU(2)-invariance (7l[Tl|. 

Compared to the previous possibility, this scheme seems to localize the degrees of freedom at the vertices of the 
graph. 

These two paths provide interesting paramctcrizations of spinor networks, relevant to discuss their (gauge-invariant) 
dynamics and to impose further symmetries such as homogeneity or isotropy, and we will quickly overview these 
constructions in the following subsections. 

B. Reconstructing Group Elements and Spin Networks 

Let us start by identifying U(l)-invariant variables. First, one notices that the 3-vectors commute with the matching 
constraints since they are invariant under the action of phases on the spinors: 

V;>* = (zt'VK'*), {M^V; 1 } = 0, Ve,/. (10) 
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Then, in order to make the link, we reconstruct the SU(2) holonomies along the edges following @, Q- We define the 
unique group element g e £ SU(2) which maps the spinor at the source vertex to the dual spinor \zl] living at the 
target vertex: 

g e \zi) = \z% g e \z°] = -We), 9e ^ ^jjf'Jf^f ■ (H) 

V\ Z e\ Z e)( Z e\ Z e) 

The dual spinor \z] is defined through the following anti-unitary map (see in appendix for more details): 

J)- «-!»>-(■?)-»■ -(+i o 

It is clear that the group element g e defined as above is invariant under multiplication of the spinors z*' by opposed 
phases cxp(±i# e ): 

{M e ,g f } = 0,Ve,f. (12) 

Following the definition of this group element, it is straightforward to check that it maps the source vector on the 
opposite of the target vector: 

g e \zl){zl\g- 1 = \z e ][z e \ > 9e > V! = ~V e \ (13) 

where g t> denotes the action of SU(2) group elements as three-dimensional rotations acting on R 3 . 

Counting degrees of freedom, the phase space C 4£; / /\J(l) E of spinor variables after imposing the matching condition 
has dimension 8E — 2E = 2 x 3E, which matches exactly the number of U(l)-invariant variables {g e , V^} that we 
have defined. Furthermore, one can check that it is possible to write the action S^lz^} defined by (|6]) in terms of 
these new variables. Indeed, we first compute: 

(zl\d t zl) + (z e \d t z e ) = (zl\d t zl) + [z s e \d t z s e ] + [z'lg^dtgelz*] = [zl\g- l d t g e \z s e ] + d t (z s e \z s e ) . (14) 

Since g e <S SU(2), then the derivative g~ 1 dtg e decomposes on the Pauli matrices a. So its trace vanishes, Tr g~ 1 dtg e = 
0, and we can decompose the matrix | on the Pauli matrices in terms of the vector 

[zl\g- l d t g e \zl] = -V e s -^Trag-'d t g e . 
Then discarding the total derivative term 9 t (z||z|) from the action principle, we finally obtain: 



s (0) K] = / dt 



- V e s ■ Trag-'dtge + ^ p e ■ (V* + g e > V e s ) + J2 K -'J2^ 



(15) 



So solving the matching constraints, we can indeed re-express our (kinematical) action principle in terms of the 
holonomy-vector variables at the U(l)-invariant level. 

Parameterizing explicitly the group elements g = exp(iau ■ a) with the class angle a £ [0, 27r] and the unit vector 

u E S2, we can express the kinematical term iV e s ■ hTrag~ l dtg e in terms of the parameters (a e , u e ) by computing the 
derivative: 



9 e &t9e = i& • I (9tOt e )u e + cosa e sma e dtu e + sin a e u e /\d t u t 

where (u, dtu, u A dtu) form an orthogonal basis of M 3 since u's norm is fixed. To derive this formula, we have used 
the expression g — cos oil + i sin ait ■ a. Let us point out that this is not equal to the naive expression given by the 
derivative of the Lie algebra element, idt(au) ■ a . 

Furthermore, we can compute the Poisson brackets of the g e , V^' variables with each other and we get after some 
slightly tedious but straightforward calculations: 

{9e,9e} « 0, {(V e s >% (V/-*)^ = 2e ijk {V?>%, {V e s ,g e } = -tg e a, {V^g e } = +iag e . (16) 

where « refers to weak equality (i.e imposing the matching condition). All brackets between variables attached to 
different edges e 7^ / vanish trivially. Finally, note that {g e -,ge} amounts to the commutators of all the components 
of the group element g e with each other. 
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This means that this reproduces the usual holonomy-flux algebra on the graph V. It also means that it is legitimate 
to consider wave-functions of the group elements (fi{g e ) at the quantum level with the 3-vectors V*'* acting as the 
left and right derivative with respect to the g e 's. This leads back to the standard quantization scheme used in loop 
quantum gravity. 

However we would like to point out that such wave- functions <p(g e ) are not holomorphic functions of our spinor vari- 
ables. Indeed, the group elements g e contains one holomorphic term in the z's and one anti-holomorphic component. 
We can nevertheless define holomorphic holonomy variables: 

G e =\4)[z s e \. (17) 

The matrix Ge still transports the spinors, but is of rank one: it maps |z|] to but sends \zf) to 0. These holomorphic 
variables Ge are still clearly U(l)-invariant. They are not SU(2) group elements anymore, but still satisfy the right 
Poisson algebra with the 3-vectors V^f'* : 

{Ge,Ge} = 0, {(K S '%(K S,t )i} = 2e ijk (V e s ' t ) k , {V e S ,Ge} = -iGeZ, , Ge} = +i *Ge , (18) 

where the first commutator vanishes exactly and not only on the constrained surface. As it was implicitly implied in 
and explicitly shown in [ll[, the Hilbert space of holomorphic wave-functions <p(Ge) is isomorphic to the standard 
Hilbcrt space of spin network functionals <p(g e ), an d the isomorphism can be realized through a non-trivial kernel 
/C(<7 e , Ge) allowing back and forth between the two polarizations. 



C. From Spinors to Twisted Geometry 



Another way to parameterize the U(l)-invariant phase space is to use the twisted geometry variables developed by 
Freidel and Speziale 0, [l2[ , which are particularly relevant to building semi-classical coherent spin network states for 
spinfoam models. 

The starting point of this approach is to define a spinor z through the unique group element which maps the "origin" 
spinor = ( 3. ] onto the normalized spinor zj (z\z). Then we decompose that group clement as the product of a 



v / 

rotation along the z-axis and a rotation with axis in the {Oxy) plane: 
\z) = ^W)n(Z) e»°* |fi) = e^^W)n(Z) |fi) = ^=5 ( \ ) <Z) = j=^ 2 {z ~l) ' (19) 

where n(Z) G SU(2) defines a rotation with axis in the (Oxy) plane 1 . This means that we have parameterized the 
spinor z with four real numbers <p, (z\z) € K 2 and Z G C, so there shouldn't be any (continuous) redundancy in this 
parametrization. We can furthermore easily check that the map is one-to-one, since we can invert this definition: 

Z=*, e^=|^|, (z\z)=zl + zl 
zo \zo\ 

The complex variable Z defines the direction of the 3-vector V. Indeed, we have: 

\z)(z\ - (z\z) n(Z)\n)(n\n(Zr\ 
which is straightforward to translate into 3-vectors: 

\z)(z\ = i(\V\I + 9-S), 
|fi)(fi| = l(I + ( T3) = i(I + e z - ; ), 



1 The SU(2) group element n(Z) is a rotation e a m with rh ■ e z = and \ih\ 2 = 1. Indeed we can expand this group element as: 

m = m x + irriy, \m\ 2 = l, 



1 m 
m 1 



which matches n(Z) for Z = % tan dm. 
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Thus we simply have V/|V| = n(Z) [> e z . Finally, it is also useful to write explicitly the dual spinor \z] in terms of 
the same variables. Using that eheT 1 = h for all SU(2) group elements h and in particular for h = n(Z), we easily 



get 2 : 



- e-^n(Z)e\n) = n{Z)ee- l ' t " J3 |fi) = n(Z) e +l0<T3 with fi = |fi] = efi = ( " ) . (20) 



Applying these definitions to all the spinors zf'*, we can write all the functions over the C iE phase space in terms 
of the new variables and Z%' . For instance, we can express the holonomy variables g e in terms of these twisted 
geometry variables. The group element g e sends back to the origin ft and then maps it to \z t e \: 

g e = niZDee-^*" 3 e"^ 3 n{Z s e )~ l = n(Z\) ee^ 3 n{Z s e )- x , with £ e = + #) , (21) 

where the variables Zp* and £ e are all U(l)-invariant. 

Finally, we can write the whole action in terms of the twisted geometry variables. First, we compute the kinematical 
term for a single spinor, (z\dtz) , in terms of the new variables: 

(z\d t z) = i(z\z)dt<t> + (z\z) {n\n{Z)- 1 d t n{Z)\n) + ^d f (z\z) . 

We further compute explicitly the derivative term in the variable complex Z G C : 

(Q\n(Z)-'d t n(Zm) = 2( \ + |z|2) = i (1 1 z|a) M, with 8 = ArgZ . (22) 
Applying this formula for z = zf '* and discarding the total derivatives, we obtain: 

S^[z" e ] = / * E A - d ^ + A * ( (1 + \Z%) dtel + (L+^W) + (d0SUrC constraints )' ( 23 ) 

where A e = (z||z|) = (z'|z*) denotes the area dual to the edge e. For the details of the Poisson algebra of the twisted 
geometry variables A e ,t; e , |^|'*|, the interested reader can find a thorough analysis in Let us just point out 

that although the variables A e are SU(2)-invariant, all the remaining variables £ e , \Zl' t \,9 s e ,t have non-trivial Poisson 
brackets with the closure constraints and are thus not SU(2)-invariant. In order to get SU(2)-invariants from the 
spinor variables or the twisted geometry variables, you can use cross-ratio observables as introduced in [22j . but we 
will not go into the details of this construction in the present work. 
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2 We can actually re-absorb the e factor in the n(Z) group element itself by switching from Z to — 1/Z. After a few algebraic manipulations, 
we get: 

n (-=J = n{Z)e~ 1 e i0<T3 , with 9 = ArgZ. 

Thus, we can re-write the dual spinor as: 

\z] = -vWe-'^nf-i) l«> = -vWn(-=) e~^+ B ^ \Vt) . 
This hints that there might be a more convenient choice of complex variable Z = IZe^ instead of Z = Re lS such that: 

We can parameterize both variables using a "boost" parameter, thus writing R = sinhr; and 1Z = tanhr;. Finally the derivative term in 
Z simply reduces to: 

{Q.\n{Z)- 1 dtn{Z)\Q,) = ill 2 dtB = ~{Zd t Z - Zd t Z) . 

We can even go further and re-absorb the area factor A = (z\z) into the definition of the complex variable. Indeed since A is real, it 
factors out of the anti-symmetric combination ZdtZ — ZdtZ which is only sensitive to the argument of Z: 

VAZdt{y/AZ) - y/AZdt(VAZ) = A(Zd t Z - ZdtZ) . 
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D. Solving the SU(2)-invariance: U(JV) Formalism 



We have up to now reviewed the path of first implementing the U(l)-invariance on the spinor phase space. This 
has lead us to the standard holonomy-flux algebra of loop gravity and to the twisted geometry parameterization. In 
this subsection, we will give a quick overview of the other alternative of first implementing the SU(2)-invariance. This 
approach has been developed in 0, @, IT3T - fl5l Il9j and been referred to as the U(iV) formalism for SU(2) intertwiners. 
Here, we will not review this whole framework, but will focus on the definition of the classical SU(2)-invariant 
observables and how to write our action at the kinematical level in terms of them. 

Focusing on the closure constraints at vertices, a natural set of SU(2)-invariant observables at each vertex v is given 
by: 



JJIV 



f 



WW), 



(24) 



Calling N v the valency of the vertex v (number of edges attached to v), the N v x N v matrix E v is Hermitian while the 
matrix F v is anti-symmetric and holomorphic in the spinor variables. It is clear that these scalar product between 
spinors are invariant under SU(2) transformations at the vertex v. Moreover the Poisson algebra of these observables 
closes: 



} 



5ehE V gf ) 



{El f ,F; h } = -i(8 eh F] g -8 eg F] h ), {E v ef ,FZ h } = -i(6 fg F* h 

{Fef> F gh) = ~i {SegE V hf - 8 eh E V gf - SfgE^ + 6 fh E^ ge ) , 



(25) 



dfhFeg) . 



Actually, the algebra of the E'-observables alone closes and the corresponding Poisson brackets actually define a u(N v ) 
algebra, thus the name "U(iV)-formalism" . One can go further with the U(JV) idea and write the E and F observables 
in terms of a single unitary matrix U v G \J(N V ): 



E v 



X v U V A*U V , 



A = 



(26) 



0jv„ 



F° = X v U v A £ t U v , 




( i 1 










0iV„-2 



where X v G M + . It is direct to write the spinors in terms of the couple of variables (A„, U v ): 



(27) 



The parameter X v is a global scale factor and measures in fact the total area around the vertex v. The advantage of 
this formulation is that the closure constraints on the spinors at the vertex v gets simply encoded in the unitarity of 
the matrix U v , as shown in 0: 



C v 



E 



i 



K){zt\-^{zt\zl) = Q & (U'^U 



I. 



(28) 



One can check that we have the right number of degrees of freedom. Around a fixed vertex v, we had started 
with N v spinors satisfying the closure constraints and up to SU(2)-transformations, which makes 4N V — 6 variables. 
Working with the unitary matrix, we must first notice that the definition of the matrices E and F are invariant under 
\]{N V — 2) x SU(2) (which is the stabilizer group of both A and A e matrices). Thus imposing this symmetry and not 
forgetting the degree of freedom encoded in A„, we have 1 + N% — (N v — 2) 2 — 3 = 4N V — 6. 

One can finally write the action defining the kinematics on the graph T in terms of these V(N) variables: 



SM[zZ] = Jdt E +iX v TrU^d t U- 1 + E Tre v ((U'^U v - I) + E Pe(E s ee - K e ), 



(29) 



where the new Lagrange multiplier Q v imposes the unitarity of U v and the matching constraints generate multiplica- 
tions by phases on the matrix elements of the unitary U v . 
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Here, we have replaced the closure constraints on the spinors by the unitarity constraints on the matrices U v . We 
can actually go further and formulate everything entirely in terms of SU(2)-invariants: 

S^[zl] = j dtJ2^r(F^d t F v +J2pe(E s ee - 2&), (30) 
where both the E v and X v should be considered as functions of the F-variables: 

A, = i5>M= ^Tr(F")tF«, 
E v = = ^-(F V )*F V , 

where the equation giving the iJ-matrix in terms of the i^-matrix is part of a system of quadratic equations relating 
the E v e} and F v ef observables 0. This means that we write our action principle entirely in the -F 11 -variables, which 
are SU(2)-invariants. However, if we consider the F^j as our basic variables, the drawback is that they satisfy specific 
constraints. Of course, one must not forget that the matrices F v are anti-symmetric, but furthermore they satisfy 
the Plucker relations (e.g. fA H, Ii~5| ) : 

ipv tpv rpv rpv , rpv rpv 

^ij^kl — t il t kj + *ik*jl ■ 

Nevertheless, keeping this in mind, the F v -variables can be very useful, since they were shown to be the natural 
variables when quantizing the SU(2)-invariant phase space in order to recover the Hilbert space of SU(2) intertwiners 
and to construct coherent intertwiner states [3, LLJ, LL5| . From this perspective, the Plucker relations define the basic 
recoupling relations for the spin-i representation. 

The insights that one has to keep in mind from this formalism are: 

• At each vertex v, the unitary matrix U v up to U(N V — 2) x SU(2) transformations encodes the SU(2)-invariant 
information about the spinors, i.e the shape of the dual polyhedron at the classical level and the intertwiner 
state at the quantum level. 

• When quantizing, we can replace the wave- functions f(z^) by wave- functions Lp{U v ) or actually by wave- functions 
f(F v ) which already implement the SU(2Vgauge invariance at every vertex. The equivalence between these 
formulations have already been studied in j7l. Il3l. [l4| . 

• There is a natural \J(N V ) action on the spinors at each vertex. Dropping the index v for simplicity's sake, it 
reads for a transformation V £ \J(N): 

U^U = VU, z e ->I e = V ef z f . 

The closure constraints commute with this U(N) action. These U(iV) transformations deform the shape of the 
dual polyhedron (or cquivalently the intertwiner at the quantum level) while keeping the total boundary area 
unchanged. And one can define coherent intertwiners which transforms covariantly under such transformations 

mill- 



E. Defining LQG Dynamics on a Fixed Graph 

Up to now, we have described the loop gravity kinematics on a fixed graph T, either in terms of spinors or U(l)- 
invariants or SU(2)-invariants. We have discussed the physical interpretation of these variables as defining a discrete 
geometry and we have introduced an action principle [z^] encoding the kinematical phase space structure together 
with the closure and matching constraints generating the U(l) and SU(2) gauge symmetries. The next step is to endow 
these kinematical structures with dynamics. The natural way to do so is to add a Hamiltonian term of our action 
and define: 

S[zl] = S^[zl]- j dtH, (31) 

where the Hamiltonian H will be a gauge invariant functional of the spinors (and possibly depending on the time 
t). This is the nice advantage of the spinorial formalism: the kinematics and phase space structure, and thus the 
dynamics, can be simply written in an action principle. 
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The dynamics defined by H will either be a Hamiltonian constraint or the evolution with respect to an external 
or gauge-fixing time parameter. We can construct H by different ways. We could discretize general relativity's 
Hamiltonian constraint on a fixed graph (e.g. [HI)' we could truncate the loop quantum gravity dynamics to T (e.g. 
[HQ), we could also extract an effective Hamiltonian evolution from spinfoam transition amplitudes (e.g. Q), but we 
can also consider all the possibilities of gauge-invariant dynamics (compatible with some notion of homogeneity and 
isotropy of the chosen graph T) and investigate which of these dynamics can be interpreted as the classical evolution 
(plus potential quantum corrections) in general relativity of certain 3-metrics. 

From this perspective, it is very important to understand what kind of 3-metrics and 4-mctrics can be generated 
by restricting oneself to a fixed graph T in loop gravity. Our point of view is that working on a fixed graph does not 
necessarily need to be interpreted as describing the evolution of discrete geometries. Indeed, one can view the fixed 
graph and the data living on it as a triangulation of a continuous geometry, a sampling from which we can reconstruct 
the whole geometry if the triangulation is refined enough and the geometry is smooth enough to interpolate between the 
discrete sampling that the triangulation has defined. The same perspective has been developed recently independently 
in [IJ], where they describe the class of continuum space geometries compatible with the finite-dimensional phase 
space structure of loop quantum gravity on a fixed graph. In other words, we consider working on a fixed graph as 
defining a mini-superspace model of (loop) quantum gravity 4 , or possibly a midi-superspace model in the case where 
we consider a certain class of graphs (e.g. 2- vertex graphs as in Q) and look at the limit with infinite number of edges. 
In this context, families of geometries defined by a finite number of parameters, as in mini- and midi- superspace 
models, can be effectively represented by the data living on a fixed finite graph, even if these geometries are not 
intrinsically discrete and even if they are not compact. For example, as we will see in the next sections, following the 
original ideas developed in [!, 0, [f| , 2- vertex graphs are well-suited to describe homogeneous and isotropic cosmology 
of the FKW type (and possibly also some of Bianchi's models). Then, having fixed the graph T defining the geometry 
on spatial slices and having identified the geometries that we wish to describe, we would like to define a proper 
dynamics of the spinor variables living on T in order to reproduce the correct dynamics of general relativity for these 
geometries, with possibly effective quantum gravity corrections. 

Thus, the interpretation that we propose of truncating the loop quantum gravity dynamics to a fixed graph definitely 
lies within mini-superspace models. We do not insist that the classical spinor data or the spin network state living 
on this graph defines the entire space itself and its geometry as in standard loop quantum gravity. But we mean that 
the kinematics and dynamics of some restricted class of geometries can be effectively modeled by the dynamics of 
classical spinor networks on some fixed graph. 



FIG. 1: The flower graph with a single vertex and L petals, on the left, and the 2-vertex graph with L + 1 edges between its 
two vertices: both define the Hilbert space of SU(2)-invariant spin network states L 2 (SU(2) i /AdSU(2)). 

A last remark that we would like to do in this section is on the role and relevance of the chosen graph T. On the 
mathematical level, the phase space of gauge-invariant observables on a fixed graph T, defined as C 4E //(U(1) E x 
SU(2) y ), turns out not to depend on the particular combinatorics of the graph T but only on the number of edges and 
vertices. To understand this, let us think in terms of the standard cylindrical functionals of loop quantum gravity, 
that is the Hilbert space L 2 (SU(2) E /SU(2) V ) of SU(2)-invariant functions tp({g e }). This space is isomorphic to 
L 2 (SU(2) B - y ) (assuming that E > V, else it is isomorphic to L 2 (SU(2)/AdSU(2)) when E = V), which only depends 



4 By mini-superspace models, we mean here a truncation of the gravity phase space to a finite number of (kinematical) degrees of freedom. 
In the context of cosmology, this is usually achieved by an appropriate symmetry reduction to a homogeneous sector of general relativity. 
But it can be more generally understood as the truncation of general relativity to a certain metric ansatz defined by a finite number of 
parameters, which then become the degrees of freedom of the mini-superspace model. That's what is achieved by working on a fixed 
graph in our context. Nevertheless, once having fixed the background graph, one can further work out some symmetry reduction, like 
the restriction to the homogeneous and isotropic sector for the 2-vertex model. 
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on the number of (independent) loops of the graph L = E — V + 1. Going further, the space L 2 (SU(2)- E /SU(2) l/ ) is 
actually isomorphic to any space L 2 (SU(2) B+TI /SU(2) y+ ™) with n > —V 1 n € Z. This isomorphism can be realized 
exactly through gauge fixing (or "unfixing") of the SU(2)-invariance at the vertices of the graph T (see e.g. [If| for 
a rigorous approach to gauge- fixing spin network functionals) . From this point of view, the Hilbert space of a fixed 
graph is always isomorphic to the space of states on a flower graph L 2 (S\J (2) L /Ad S\J (2)) or on a 2-vertex graph 
L 2 (SU(2) L+1 /SU(2) 2 ). Thus we can always write, both at the classical and quantum levels, any gauge-invariant 
dynamics defined on T on a corresponding flower or 2-vertex graph. 

The natural question is then why should we bother about using more complicated graphs involving more vertices 
and a more complex combinatorial structure. Our point of view is that the combinatorial structure of the graph T 
provides us with an implicit vision of the space geometry. For instance, through the implicit notion that a vertex 
of the graph represents a physical point of space (or region of space) and that edges define directions around these 
points, the graph's structure will define our notions of homogeneity and isotropy. Therefore, even though we can 
translate any gauge-invariant observable or Hamiltonian operator from our potentially-complicated graph T to any 
graph with the same number of loops, the concepts of locality, homogeneity and isotropy of our dynamics crucially 
depend on our original choice of T. From this viewpoint, different graphs will admit a different implementation of the 
symmetry reduction (to e.g. the homogeneous sector), that in turn will lead to different mini-supcrspace models. 



Now that we have introduced and reviewed the spinorial formulation of the loop gravity phase space on an arbitrary 
fixed graph T, and discussed how to recover the standard observables -holonomies, fluxes, twisted geometry variables, 
U(iV) observables- in this framework, we will now specialize to the case of the 2-vertex graph (see figH]). Indeed, there 
is serious evidence that truncating the loop quantum gravity dynamics to this very simple graph and implementing an 
appropriate symmetry reduction leads to Fricdmann- Robertson- Walker (FRW) homogeneous and isotropic cosmologies 
|3j-l7[. In this section, we will study the kinematics on the 2-vertex graph, clear up in this simpler context the 
geometrical meaning of the variables defined in the previous section and explain how to reduce to the homogeneous 
and isotropic sector following the ideas developed in @, 0| • We will tackle the issue of the dynamics and how the 
Hamiltonian generates the four-dimensional FRW metric in the next section. 



Let a and j3 denote the two vertices of the graph as on figlU Following Sec. HI we denote by and Wi, with 
i = 1,...,N, the collection of spinors living on the N edges respectively attached to the vertex a and to the vertex fj. 
These spinors are subject to the closure constraints, 



II. 



CLASSICAL PHASE SPACE ON THE 2-VERTEX GRAPH 




FIG. 2: The 2-vertex graph with N edges linking the two vertices a and /3. 



A. Kinematics on the 2- Vertex Graph 




(32) 




with X(z) = | J2i( z i\ z i) = h Si r i z i)\i an d to the matching constraints, 



(zi\zi) = (wi\wi) Vi = l.JV <=>• 1^)1 = 1^)1. (33) 
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In particular, the total boundary area is of course the same seen from the two vertices, A = X(z) = X(w),. 

The action principle summarizing the phase space structure can be cquivalcntly written as a functional of the 
spinors Zk,Wk or of the vectors V k S K 3 and group elements g k & SU(2) : 

N 

S {0} [z k ,w k ] = / dt^2-i((z k \d t z k ) + (w k \d t w k )) + (z k \A a \z k ) + (w k \Ap\w k ) + pk((z k \z k ) - (w k \w k )) (34) 

. JV N N 

= / dtY, V(z k ) ■ l -Trag^d t g k + A a ■ + Ap • ^ V{w k ) + p k ■ (V(w k ) + g k > V{z k )). 

J k=l k=l k=l 

We can further re-express this action in terms of the twisted geometry variables or in terms of the SU(2)-invariant 
E a 'P and F a ^ living at both vertices. 

The interpretation of the 2-vertex graph as a discrete geometry is clearly of two polyhedra glued along their N 
faces with matching areas. Nevertheless, the two polyhedra can still have different shapes and this can be interpreted 
as local degrees of freedom living at each vertex. We will see below that constraining the two polyhedra to have the 
same shape naturally projects us onto the homogeneous and isotropic sector of the phase space. 



B. Symmetry Reduction to a Isotropic Cosmological Setting 

In @, 0], a symmetry reduction of the 2-vertex graph model was proposed by requiring a global U(iV) invariance of 
the spinor networks. Indeed, this invariance leads to an homogeneous and isotropic model parameterized by just one 
degree of freedom. Let us see how this reduction works. It is best understood when written in the U(iV) observables 
which are SU(2)-invariant. Reminding the definitions given above in section ITPl we have two sets of SU(2)-invariant 
quantities defined at the vertices of the graph: 

Ef j = {z i \z j ), F% = [ Zi \ Zj ), E^( Wl \w 3 ), F| = KK). 

The matching constraints on the N edges are simply expressed in terms of these observables: 

Vi, Mi = E%-E& = 0, 

and they generate an invariance under XJ(1) N as said in the previous section. The point noticed in [(| 0] is that this 
symmetry can be naturally enlarged from 11(1)^ to U(7V). Let us indeed introduce: 

Vi,j, £ ij = E§-E p ji , £u=Mi. (35) 

Considering the Poisson bracket of the £?-observables, it is straightforward to check that these new SU(2)-invariants 
form a U(AT) algebra: 

{£ij,£kl} = —i (fijAi - SuSkj) ■ (36) 
Let us further impose these constraints on the phase space of our spinor variables living on the 2-vertex graph: 

f N 

S h}mo[ z ki w k\ = I dt^2-i({zk\d t Zk) + (w k \d t w k )) + (z k \A a \z k ) + (w k \Ap\w k ) + p kl £ k i, (37) 

^ k=l 

where the matching constraints M. k are taken into account as the diagonal components (when k = I) of the U(iV)- 
constraints £ k \. 

This new set of constraints £ geometrically means that the polyedra dual to the two vertices a and /3 are identical. 
Indeed requiring Ef^ = E^ for all couple of (possibly identical) edges i,j implies (by taking the norm squared) that 

all the scalar products between 3- vectors Vi-Vj at the two vertices are equal. This tells us that these U(iV) constraints 
impose homogeneity: the intertwiner states at the two vertices are identical. 

Coming back to the spinors, the constraints read (zi\zj) = [wi\wj] for all i, j's, which implies after a little algebra 5 
that the spinors \zi) are equal to the spinors \wi) up to global U(2) transformation (it's an equality on scalar products) 
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(see e.g. @, [HI). Since both sets of spinors z% and Wi are anyway defined up to SU(2) transformations generated by 
the closure constraints at the vertices a and /3, this means that the spinors are exactly equal to the spinors \wi] 
up to a global phase: 

Vfc.Z, £ki=0 => Vi \w t ] = e* \zi) . (38) 

This means that the unitary matrices U a and U 13 are also equal up to a phase, U 13 = e l ^U a (up to SU(2) x \J(N — 2) 
transformations). We can also translate this condition in the SU(2) group elements g k living on the edges. Just being 
careful that e 1 ^! ^ SU(2), the holonomics read: 

= e^\z k )(z k \+e^\z k ][z k \ 
9k (z k \z k ) 

thus being the diagonal SU(2) matrix (e ! ^,e -1 ^) in the \z k ), \z k ] orthonormal basis. 

Finally, putting back the solution \ wi] = e 1 ^ \z%) in the action (|53|) with U(iV) constraints, the action greatly reduces 
and we simply get: 

S ( hLo[^,w k ] = -2 y"d* XM = S^[X,4>] . (39) 

The only remaining degrees of freedom are thus the total area A = 2A and its conjugate angle <p. Thus we have 
reduced the 2-vertex graph to its homogeneous and isotropic sector. 

By isotropic, wc mean that the only relevant degree of freedom is the total area and not the individual areas 
A k = (z k \z k ) = (w k \w k ). So the individual areas are not required to be the same on all edges. On the other hand, it 
means that the conjugate variable defined by the angle <f> is the same for all edges. We can go further in understanding 
the U(7V)-constraints £ k i leading to homogeneity. Indeed, they generate the following U(7V)-transformations on the 
spinors: 



(Uz) k = Uki zi, w k -> (Uw) k = Via wi . (40) 



It is this huge gauge-invariance that kills the dependence on degrees of freedom living on each edge and effectively 
reduces the action to its isotropic sector. 

We would also like to point out that the angle <j) in the symmetry-reduced sector matches the angle £ parameterizing 
twisted geometries, which provides it with a physical interpretation as related to the extrinsic geometry, i.e. related 
to the embedding of our spatial slice -the 2-vertex graph- into the 4d space-time. 

Therefore, to conclude this subsection, once we have accomplished the symmetry reduction by the U(./V)-constraints, 
we are left with homogeneous and isotropic 2-vertex spinor networks, which is described by the holonomies e 1 ^ along 
all the edges of the graph and by the total boundary area A = 2A around both vertices. 

Our final comment is that the resulting homogeneous and isotropic sector is independent from the number of edges 
N of the original 2-vertex graph. This is normal since imposing isotropy means getting rid of the dependence on the 
individual edges. On the other hand, the full theory and anisotropy should crucially depend on the allowed number 
of the edges N. 

C. On the Choice of Complex Variables 

As we have seen, in the isotropic sector, the phase space reduces to a one-dimensional system described by the total 
area A = 2A and its conjugate angle <p: 

{A,4>} = 1. 



5 For example, we can translate the constraints on the E's into constraints on the F's using the fact that + \z\[z\ = (z\z)I for any 

spinor z: 



Vi,j,k, F ki F k j — EijE kk — E ik E k j 
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Thinking of these as conjugate variables, it is usual to choose a complex structure defined by the variable £ = 
(A - icf>)/V2: 

{(,(} = *■ (41) 

This is the standard choice done in the context of twisted geometries 0, [l2[ and in the spinfoam cosmology approach 
of Q which uses the twisted geometry phase space to describe its boundary states. Another choice, better suited to 
the spinor approach used here and to the coherent spin network states based on the spinorial framework fioj , is the 
complex variable z — \J2A e~ l i : 

{z,z} = i. (42) 

This is also a canonical choice of complex structure on this two-dimensional phase space. This is actually the choice 
that we will make when analyzing the dynamics of the 2-vertex graph phase space induced by the spinfoam amplitudes, 
as we will see in section HVl 

The advantage of z over f is that it reflects the original definition of the variables. Indeed, A is defined as a real 
positive number and <f> is defined as an angle (modulo 2n) from the holonomy. Then z = \JlA e~ l 2 is a more natural 
complex number in order to define the unitary matrices U a '^ . On the other hand, using the complex structure defined 
by £ = (A — i0)/v2, and quantizing the system in terms of this variable, means that one has to restrict by hand A to 
be positive and <p to be bounded. In simpler terms, z £ C parameterizes exactly our phase space taking into account 
the constraints A 6 R + , <fi € [0, 2ir] while ( £ C leads to some redundancies that have to be dealt with in a nontrivial 
way at the quantum level. Finally, as we will see in section HV[ the spinfoam transition amplitudes in the isotropic 
and homogeneous sector will naturally be holomorphic functions of z. 



D. Beyond the Isotropic Sector 

In the present work, we focus on the study of the homogeneous and isotropic sector and its dynamics both at 
the classical level through the definition of classical Hamiltonians as in the following section IIIII and at the quantum 
level through the analysis of the spinfoam transition amplitudes in section HVl Nevertheless, from the perspective of 
cosmology, it is necessary to go beyond this simplistic model and study the departures from isotropy and homogeneity. 
In particular, we should study the dynamics of the inhomogeneities and their feed-back on the global homogeneous 
sector in order to compare to current measurements in cosmology (on the cosmic microwave background for instance) . 
Extending the analysis of our 2-vertex model to the whole Hilbcrt space, beyond the U(iV)-invariant sector, would 
thus allow to test its relevance and validity as a quantum mini-superspace model for cosmology. 

One goal would be to parameterize efficiently the reduced gauge-invariant phase space on the 2-vertex graph and 
understand how to project the phase space variables onto the spherical harmonics and define a multi-pole expansion 
of the corresponding degrees of freedom. For instance, the number of edges iV is irrelevant to the kinematics and 
dynamics of the homogeneous and isotropic sector since requiring the U(iV)-invariancc kills any dependence on N. On 
the other hand, N will play a non-trivial role out of the U(./V)-invariant sector as a cut-off in the number of degrees 
of freedom. It would then be interesting to understand how big should N be to model our observed cosmology and, 
for example, allow anisotropy as seen in the cosmic microwave background. 

Actually, an analysis of the 2-vertex model beyond the isotropic sector has been carried out for the case N = 4 
p6j . In that work, the 2-vertex model is regarded as a triangulation of the whole space, that is thus a three sphere. 
The reduced gauge-invariant phase space, which in this N = 4 case is formed by twelve degrees of freedom, has 
been then identified with the Bianchi IX model (homogeneous and isotropic model with the spatial topology of S* 3 ) 
plus perturbations. Since the Bianchi IX model is described by six degrees of freedom, those perturbations account 
for the remaining six degrees of freedom. Using the fact that S 3 is isomorphic to SU(2), in |26| a expansion of the 
perturbations in terms of Wigner matrices is considered, and those six remaining degrees of freedom are identified 
with the six components forming the diagonal part of the lowest integer mode in that expansion, though the authors 
do not give a justification for such identification. 

We will not presently investigate in detail how to parameterize the full phase space out of the U(iV)-invariant sector 
for general N but merely discuss the possibilities from our perspective. We postpone the full analysis to future work. 
The degrees of freedom of the 2-vertex graph can be understood from the point of view of both the edges, as the 
holonomies defining the curvature between the two vertices, or both the vertices, as the internal geometries of the two 
vertices and the correlations between them. If one decides to focus on the holonomies, it seems natural to use the 
holomorphic and anti-holomorphic components of the holonomies around the loops of the graph (expressed in terms 
of the E and F variables) as introduced in Q and to look for a way to combine them in order to get a finite Poisson 
algebra encoding the whole reduced phase space. The alternative is to start with the unitary matrices U a '^ describing 
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the geometry of each vertex. Let us look a bit more into how to construct gauge invariant SU(2) observables from 
these matrices. 

First of all, as seen in section llD| the observables must be invariant under the right action by SU(2) x U(iV — 2) 
at both vertices: 

U a -^U a V a , U f3 ^U f3 V f3 , V a , V 13 € SU(2) x U(A — 2) . 

It is thus natural to make U a and U 13 act on SU(2) x U(A — 2) vectors. Following the previous work [l3|, [H|, we 
choose the irreducible representations of U(iV) whose highest vector is invariant under SU(2) x \](N — 2). These are 
labeled by an integer j£ N and correspond to Young tableaux with two horizontal lines of equal length J. Moreover 
they have been shown to be exactly the Hilbert spaces of SU(2)-intertwiners (at fixed total area J) resulting from the 
quantization of the spinorial phase space 0, EH, EIJ- Let us ca ^ l-^ ^n) the highest weight vector of these irreducible 
representations. As shown in (T3], acting with the matrix U a € U(iV) on \J, f2jy) generates a coherent intertwiner 
state peaked on the classical phase space point with the spinors given by the first columns of the unitary matrix U a 
according to the classical definition ((27)) : 

U a \J,fl N ) = \J,{z?} i=1 .. N ), (43) 

where the zf satisfy the global normalization condition J2i( z ?\ z i) = 1- The highest weight vector can be identified 
as the bivalent case, where all spinors vanish except the first two: 

''■«">-''.{(;).(;).(::).■•}> <-» 

Then acting with the unitary matrix U a on this state, we generate all coherent intertwiners. We similarly define the 
lowest weight vector by acting with the duality map on the spinors: 

\J,&n) - 

More details on this construction can be found in [lj, [15( . Then the scalar product between the states in a and ft is 
invariant under V(N) (acting as U a -> UU a , U p -> UU 13 ) and can be computed exactly: 

(?f2jv, J\ t U^U a \Slif, J) = <4{«f}IW}> - ^det£|*f)[*fll . (45) 




The first remark is that the zf'^ are just our spinors Zj and Wi up to the normalization factor yf\ defining the total 
area. Then the expression simplifies greatly assuming that the Zj's and w, 's satisfy the closure constraints and that 
we are in the isotropic sector with \wi] = e l ^\zi) and we simply get: 

{n N ,J\ l U p U a \n Nl J) = (e-^y. (46) 

This scalar product can thus be considered as a definition of the conjugate angle </> in the full phase space. 

The second remark is that the scalar product (f2jvj J\ t U^U a \Q,N, J) is exactly the evaluation at the identity of the 
2-vertex spin network state with coherent intertwiners attached to the vertices a and ft. It is then straightforward to 
obtain an overcomplete basis of gauge-invariant observables by considering the evaluation of this spin network state 
on arbitrary group elements gi € SU(2) living on every edge of the graph. This amounts to inserting these group 
elements in the scalar product expression between the two U(iV) matrices. Such evaluations are by definition both 
SU(2)-invariant at the vertices and U(l)-invariant along the edges. These evaluations and their Poisson algebra might 
not be easy to compute, but this procedure still hints towards a parametrization of the gauge-invariant phase space 
in terms of U(iV) representations. From this viewpoint, it might be interesting to investigate the relation between the 
Peter- Weyl expansion in U(iV) representation and a multi-pole expansion of the observables in our phase space. We 
would then expect to recover the continuum limit (i.e no cut-off in j in the spherical harmonics expansion) as N is 
sent to infinity. 

A detailed analysis of the full phase space is postponed to future investigation of anisotropy in the 2-vertex model 
for cosmology. 
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III. CLASSICAL DYNAMICS ON THE 2- VERTEX GRAPH 



As said before, the above kincmatical setting, namely the U(N)-invariant spinor network defined on the 2-vertex 
graph, seems suitable to model effective dynamics for Friedmann-Robertson- Walker (FRW) like cosmologies. In this 
section we will face this issue from the canonical point of view, by adding an appropriate Hamiltonian to the action, 
as it was done in Q. We will analyze whether the resulting model can be understood as an effective FRW model, 
that introduces corrections coming from the discrete theory. 



A. Hamiltonian for FRW Cosmology 

Before describing the effective dynamics that arises naturally from the two-vertex model, let us first review the 
dynamics of FRW models as dictated by general relativity, so that we can make the link between both theories. 
The FRW model represents homogeneous and isotropic spacetimes. Their 4- metric is given by the general form 

ds 2 = -N 2 {t)dt 2 + a 2 {t)dS 2 , dS 2 =°q ab dx a dx b . (47) 

Here Af is the lapse function and a(t) is the scale factor, that we consider with dimensions of length so that the 
coordinates x a are dimcnsionless. The indices a and b go from 1 to 3 and denote spatial indices. Finally, °q a b stands 
for the 3-metric of a 3-dimcnsional space of uniform curvature. There exist three such spaces: Euclidean space, 
spherical space, or hyperbolic space. In (dimcnsionless) polar coordinates (r, 9, <j)) we simply have 

r 2 

dS 2 = — + r 2 (d9 2 + sin 2 6>d0 2 ), (48) 

1 — kr z 

where k is a constant representing the curvature of the space: k = 1 for the sphere, k = — 1 for the hypcrboloid, and 
k = for the Euclidean space. 

In order to get the dynamics of a generic FRW model, the easiest way to proceed is to take the Einstein-Hilbcrt 
action for general relativity with cosmological constant A, 

1 



Sgr = J dt _/ s d 3 xy/d^(R[g] + A), (49) 

and particularize it to the FRW metric given in Eqs. (|47 |) -(j48 |) . After employing a 3+1 ADM decomposition to express 
it in the Hamiltonian formalism, the result for Sgr[9frw] is 

Sfrw = I dt [ d 3 xV^[n a d t a-N-C] : C = J^LA _ + 3. (50 ) 

J 7S 6 a 07T(_r 07r(_r 

In the above expressions G denotes the Newton constant, °q denotes the determinant of the fiducial metric °q a b, and £ 
stands for a 3-dimcnsional spatial slice of the space-time. In presence of homogeneity, since all the functions entering 
the action are homogeneous, the spatial integral over the slices S diverges for the non-compact spatial topologies 
corresponding to k = 0, — 1. Nonetheless, one usually avoids that divergence by restricting the above spatial integral 
just to a finite cell V with finite fiducial volume 



V = / d 3 x^q~. (51) 



V 

Indeed, because of the homogeneity, the dynamics of the whole spacetime is the same in every point and then a finite 
region is enough to capture the dynamics. 

Note that the lapse Af is a Lagrange multiplier which imposes the Hamiltonian constraint C — 0. Actually, since 
this model is homogeneous, the spatial diffcomorphism symmetry of general relativity trivializes and there is only left 
the invariance under time reparameterizations, which is implemented by the Hamiltonian constraint. 

Remarkably, the FRW cosmologies, being a solution of the equations of motion obtained by varying the reduced 
action (f5T)]) . are also a particular solution of the Einstein equations (that are obtained by varying the full action (|4"9")1). 
In other words, the FRW model is a symmetry reduction of general relativity, more concretely, its homogeneous and 
isotropic sector. 

In the vacuum case under study it is very easy to obtain the trajectories in phase space, 7r a (a), using the Hamiltonian 
constraint C = 0. They are given by 



tt q (a) ± — - \/3Aa 2 - 9k. (52) 
47rG 
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Vanishing cosmological constant Negative cosmological constant Positive cosmological constant 




| » k=0 k=-l| k=-l | — k-0 k-1 k=-l| 



FIG. 3: Trajectories n a {a) in the different vacuum FRW models. In these plots we have used the conventions 4irG — 1, and 
A = ±1 for the cases with non-vanishing cosmological constant. 



In fig. [3] we show them just for positive ir a . 

This vacuum case is not very interesting as long as there are no degrees of freedom. In fact, the geometry is totally 
fixed by the Hamiltonian constraint C = 0. As a consequence many trajectories are rather boring, like e.g that of the 
fiat case (k = 0) with vanishing cosmological constant. Other trajectories are not even real, as happens for the cases 
{A = 0, k = 1}, {A < 0, k = 0}, {A < 0, k = 1}. In order to have more interesting phenomenology one should add 
matter to the model. Actually, already the simplest form of matter, a minimally coupled massless scalar, provides a 
quite non-trivial evolution. However in this paper we will just focus on the vacuum case, as a step to start with in 
our derivation of effective dynamics for cosmology using the 2- vertex graph. In the next section, we will derive this 
effective dynamics and discuss if it can be mapped in some regime to the cases occurring in the FRW model. We 
leave for the future the study of non-vacuum effective models. 



B. Cosmological Hamiltonian on the 2- Vertex Graph 

Let us now analyze the addition of an appropriate Hamiltonian to the kinematical action Q39[) . This Hamiltonian 
must be, first <S'J7(2)-invarint, so that the gauge invariance of the theory is preserved, and second £/(iV)-invariant, 
so that it is homogeneous and isotropic. In this way, the resulting dynamics will be consistent with the kinematical 
setting, and even more, it may be regarded as generating the reduced (homogeneous and isotropic) sector of the full 
theory 6 . In order to construct such an ansatz for the Hamiltonian, the simplest SU(2) invariants to consider are E a 
and . The next simplest possibility is to consider 5C/(2)-invariant observablcs acting on elementary loops made 
of two edges. These observables are of the form E^E^f, F"^F^, and F^F^. Now, the symmetry reduction to the 
homogeneous and isotropic sector implemented by the U(N )-invariance, reduces the above 517(2) invariants to the 
[/(AO-invariant terms TrE a = TrE 13 , Tr t E a E f) oc (TrE a ) 2 , TrF a F 13 , and TrF a F^. Made up of these terms, we 
will consider the following ansatz for an action with non-trivial dynamics, on the U {N )-invariant two- vertex spinor 
network 

Shorno[zk, Wk] = S^ mQ [z k , W k ] - J dtff (H homo [z k , W k ] - H Q ), 

H homo [z k , w k ] = 1 °Tr t E a E f} + 1 +TrF a F 13 + y^TrF a F^ + ^[Tr*£ Q ] 3 (53) 

Here 7 , 7 + , 7", and 7 1 are some real coupling constants. On the other hand M is a Lagrange multiplier imposing the 
Hamiltonian constraint Hhomo[z k , w k ] — H Q = 0, with H Q a real constant accounting for the fact that in the quantum 
theory the energy of the fundamental state could be nonzero 7 . In the following we will take N = 1. The above ansatz 



6 As the FRW Hamiltonian does for general relativity. 

7 In the same manner as the energy of the fundamental state of the harmonic oscillator or that of the hydrogen atom is not null. 
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is actually that studied in Q plus an additional term, proportional to [T^E ] 3 , that it is related with a cosmological 
constant term, as we will see in a moment. 

If we remember the form of the matrices E v and F v in terms of and U v , and that the f/(7V)-invariance implies 
jjp — e i<t>ua anc j \ a — Xp = X, the above action reduces to the homogeneous and isotropic action 



2/„,o „,+ „-2i4 



7 - e 2# + 7 i A ) (54) 



S[X, <j>] = -2J dt\d t 4> - J dt(H[X, <p] - Ho), H(X, 0) = 2A 2 ( 7 ° - 7 

We will further choose 7+ = 7~ = 7/2 so that the Hamiltonian is real and given by 

H(X,<P) = 2A 2 [ 7 ° - 7 cos(20) +j 1 X}. 

We can now compare S[X, cj)] with the action of FRW given in Eq. ((50)) . In order to do that, let us express Sfrw 
in our variables (A,0). Since A has the interpretation of an area, we can take the following canonical transformation 
to relate the variables (a, 7r a ), with Poisson bracket {a,7r a } = 1/V , with our variables (A,</>), with Poisson bracket 
{A,0} = l/2: 



(55) 



V, 



1/3 



IVo 



2/3 



VXcj>, 



(56) 



where I serves as a unity of length defining the relation between our dimensionless area A with the dimensionless scale 
factor a/7. As before, V a represents the (dimensionless) volume of a finite cell V, measured with respect to the fiducial 
metric °q a b of the 3-dimcnsional spaces with constant curvature. Namely, we regard the graph as dual (not to the 
whole space but just) to a finite region V, which it is enough for determining the dynamics of the space because of 
the homogeneity. 

Employing the above canonical transformation we obtain the Hamiltonian of the FRW model in our variables. It 
has the form 



FRW 



= V n NC = 



32ttG 



x 2 4? 



3KI 



;A 2 + 



A/ 3 



A ! 



8ttG 8ttG ' 



(57) 



where, for convenience, we have chosen N = A 3 / 2 , and we have defined K = kV 2 ^ 3 '■ Furthermore, the trajectories in 
phase space in terms of these variables are now given by 



0(A) = ±- 



r 2 



V3A/ 2 A- 9A". 



"4(4ttG) 

The positive and real trajectories (shown before in fig. [3]) have now the graphic shown in fig. |U 



(58) 



Vanishing cosmological constant 



Negative cosmological constant 



Positive cosmological constant 



k-0 — — k--l| 
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FIG. 4: Trajectories 0(A) in the different vacuum FRW models. In these plots we have used the conventions V = 1, 4-7rG = 1, 
and A = ±1 for the cases with non- vanishing cosmological constant. 



Now in the limit (f> — > we can do the approximation cos(2</>) k1- 2<p 2 and identify our Hamiltonian H with that 
of the FRW model. Indeed, upon that approximation both Hamiltonians agree by doing the identifications 

8ttG 3KI 8nG , A^ 3 



7 



3/ 3 



16ttG' 



7 



3Z 3 ' 



7 



16ttG' 



(59) 



20 



As we said before, the term with coupling constant 7 1 represents a cosmological constant term. On the other hand, 
the other two terms, with coupling constants 7 and 7 account for the curvature term and for the other term, the 
only term remaining when there is neither cosmological constant nor curvature. 

Outside the regime in which the above approximation is valid, the Hamiltonian of the two-vertex model with the 
above identifications reads 



H(\,< 



2\ 2 



IttG ZKl 



3/ 3 16ttG 



8ttG ,„ jS A/ 3 ' 
■cos(20) + — — —A 



3/ 3 



16ttG 



(60) 



Alternatively, we can express it in the variables {a, 7r a } commonly employed in cosmology, in which case is given by 



H{a,* a ) = 2V^ 



8nG 3KI \ 8nG 

! H 

3/3 T 16ttG J 3/3 



't/ 1 / 3 ;2 1 
Vo t 7T 



A * 3 F 2/3" 2 

16ttG I 2 



(61) 



We would like to emphasize that here we are pushing forward the identifications of Eq. (|59|) further from the regime 
in which <j) — ¥ 0, where they have been obtained. 

Our goal is to analyze whether the Hamiltonian of the two-vertex model -ff (A, </>) can be regarded as an effective 
Hamiltonian for the FRW models. Upon this interpretation, this Hamiltonian introduces corrections to the results 
predicted by general relativity. On the one hand we will have corrections coming from the gauge invariance of loop 
gravity, and they are negligible whenever the angle (f> approaches a vanishing value. On the other hand we are 
considering that the Hamiltonian of the two- vertex model is fixed to be equal to some constant value H that does 
not necessarily vanish, unlike in general relativity. This fact will also induce corrections. 

Before doing the comparison between our two-vertex model and the FRW models, we first need to study the phase 
space trajectories resulting in the two- vertex model. Using the constraint H(X, </>) = H D we easily get 



4>(X) = — arccos 



/ 



4ttG 



9K 3AZ 4 



/ 3 3H 
4ttG 4A 2 



(62) 



These trajectories are drawn in fig. [5j The graphics only show the trajectory for </> € [0, tt/2]. In the range 4> £ [tt/2, tt] 
the trajectory is given by the specular image of that in the previous range with respect to the edge <f> = 7r/2, and for 
bigger range it is periodic with period 7r. We show them for different values of H . This constant is usually bounded, 
cither from below or from above, as shown in table HI We see that indeed it can not be zero for many cases. 





k = 


k = 1 


k = -1 


A = 


Ho < 


Ho < 


Ho 


A= -1 


Ho < 


Ho < 


Ho < 1.9 


A = 1 


Ho > -11 


Ho > -42 


Ho > -0.9 



TABLE I: Ranges of H a . 



Let us now compare the trajectories of the FRW model in the six cases with real solution (fig. 0]) with the analog 
cases in the two-vertex model: 

• A = 0, k = 0: The trajectory of the two-vertex model agree with that of the FRW model for increasing values 
of A. Moreover A is bounded from below. The bound depends on the particualr value of H Q and it is always 
bigger than zero. 

• A = 0, k = —1: For increasing value of A the trajectory tends to a constant value (independent of H Q ), as in 
the FRW model. This value is 

x 1 r, ( 1 \ 2 $ k 

— — arccos 1 + — — 

V 2 [ \AttG) 8 

and agrees with that of the FRW model, <fi — 3/ 2 V 1 ^ 3 /167rG, for an appropriate value of the length parameter 
/. As in the previous case, A is strictly positive and its minimum depends on the particular value of H a . 
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FIG. 5: Trajectories 0(A) in the different 2-vertex models. In these plots we have used the conventions I = 1, V a = 1, 47rG = 1, 
and A = ±1 for the cases with non- vanishing cosmological constant. 



• A < 0, k = — 1: Only for H Q > 0, (f> is bounded from above, as happens in the FRW model. For H a = the 
trajectory is quite similar to that of the FRW model, approaching it as A tends to its maximum value. More 
interestingly, for values of H a slightly bigger than zero, the trajectory still agrees with that of FRW for the 
maximum value of A, but it deviates from the FRW trajectory as A decreases, in such a way that A is bounded 
also from below and it never vanishes. 

• A > 0, k = 0, ±1: In these three cases, in the two-vertex model A is bounded both from below and from above, 
while in the FRW model A is not bounded from above. 

This analysis points out the limitations of the present two- vertex model to provide an effective cosmological model. 
Actually, it is not suitable to modelize the FRW models with positive cosmological constant, since the area of the 
cell under study can not increase arbitrarily, unlike in the FRW model. However, the other FRW models admit their 
analog in the two- vertex model. In these cases, there exists an effective two-vertex model that introduces corrections 



22 



to the FRW trajectories as the area decreases, in such a way that the area turns out to be strictly positive. The 
corrections are then unimportant in the classical regime of large areas, as desirable, since in this regime any effective 
theory should agree with general relativity. The fact that in the effective model the area has a positive bound resembles 
the results obtained in LQC, where the area (or the volume) suffers a bounce so that it never vanishes (see e.g (27l - [29j ). 
Indeed, in Q was already pointed out the possible similarity between the effective dynamics of the two- vertex model 
and that of LQC. In Appendix [B] we review the effective dynamics of LQC and formulate it in our variables to make 
clear the comparison. 

As said before, this analysis is just a starting point in the derivation of effective cosmological models from loop 
gravity formulated on a fixed graph. Already before the introduction of any kind of matter, we see that our 2- 
vertex graph model does not provide any effective model for homogeneous and isotropic cosmologies with positive 
cosmological constant. The problem does not lead in the fact that the area can display a positive bound, feature 
which in turn is a consequence of the dependence of the Hamiltonian in the variable (f> through the cosine (which is 
a bounded function), but rather in the fact that the semiclassical limit fails to be the correct one, since the large 
area limit does not necessarily corresponds to — > 0. Therefore, if we want to obtain a successful effective model for 
all possible homogeneous and isotropic cosmologies, we need to improve our approach. In the canonical framework 
different possibilities seem to be at hand: 

• We could redefine the canonical transformation ([55)1 doing 

^a^^a = 7T a + /(A) (63) 

such that the canonical commutation relation {a,7f a } = 1/V is preserved. Such a modification would change 
the phase space trajectories, but since the function /(A) is arbitrary we could try to choose it conveniently, such 
that the resulting model indeed success in providing any effective FRW model 

• So far we have assumed that our variables (A, </>) are canonically related with the common ones (a, 7r a ). However 
a deeper understanding of our model makes us think that our variables, coming from a discrete theory (that in 
principle is based on a quantum theory), may not be canonically related with the classical ones (a,7r a ), but the 
canonical relation is only approximately recovered in the regime </>—>■ 0. Furthermore, related with this issue it 
is the fact that the identifications given in Eq. (|59|) could be non- valid in the regime where the approximation 
cos(20) 1 — 2(j> 2 is not longer true. The physical meaning of the coupling constants 7 , 7, and 7 1 would be then 
different for the one assumed in our previous analysis. In consequence, this idea could allow us to drastically 
change our two- vertex model to find a successful link between it and the classical FRW cosmologies 8 . 

• Bringing to our framework the ideas of LQC, other possibility is to modify the physical meaning of the angular 
variable defined out of the holonomies of the loop formalism. This can be done by rescaling the variable by a 
function of the area, 

+ -+ + ~WY (64) 

as it is done in the improved dynamics of LQC [2^], where one chooses /(A) cx \/A, or in the lattice refinement 
approach [301 ] . where a more general rescaling is considered, but still of potential form /(A) oc A a . In this way 
the new angular variable <p is no longer canonically conjugate to the area but to some function of it, for instance 
the volume in the case of the improved dynamics of LQC. Note that after the rescaling the large area limit 
corresponds to the limit 4> — ► 0, as desired. 

• Furthermore, other possibility is to work with a different graph, which admit a different reduction to the 
homogeneous and isotropic sector. Actually, it would be very interesting to obtain a model in which the global 
variable after the reduction were the volume instead of the area, as in the improved dynamics of LQC (see also 
the discussion in Appendix [Bj) . From our experience with LQC, such a result suggests that the corresponding 
reduced model would succeed in providing a desired FRW effective model, and not only for vacuo but also for 
cosmologies coupled to matter. 



It is worth to point out that a similar strategy could be investigated in the framework of LQC. There, the initial quantum dynamics 28] 
failed in providing a successful effective dynamics for FRW models, and then it was improved |29H by modifying the physical meaning 
of the angular variable defined out of the considered holonomies. Applying our line of reasoning, one could try to make the initial LQC 
dynamics succeed by modifying the commutation relations between the basic phase space observables, such that it is canonical only in 
an approximate way in the semiclassical limit. 
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In summary, just because of the limitations of the two-vertex model, pointed out from the analysis of the phase 
space trajectories, we cannot rule out this model in our aim of modelizing effective FRW cosmologies with it. The 
model formulated in the 2-vertex graph is not as simple as it seems, in the sense that we may not yet understand 
completely the exact physical meaning of the variables (A, cj>) that describe it, as suggested by the first three points 
listed above. We will postpone the investigation of those points for future work, when we will also take into account 
the coupling with matter, in order to study true dynamical cosmological models. Instead, in the next section, we will 
rather look for the derivation of effective FRW models from loop gravity formulated on the 2-vertex graph following 
the ideas of spin foam cosmology Q. 

IV. SPINFOAM DYNAMICS 

The spinfoam framework defines a path integral formalism for quantum gravity, which allows to compute well- 
defined transition amplitudes for spin network states (see e.g. for a review). We propose here to use it to define the 
transition amplitudes between coherent spin network states peaked on the classical spinor network data and derive 
from this an effective classical dynamics for the spinor networks taking into account the spinfoam quantum gravity 
effects. 

We will not review the spinfoam framework, and we will assume that the reader is familiar with the various spinfoam 
constructions and methods. We will only introduce the necessary concepts for our derivation and refer to the known 
literature on the subject for the details. 

A. The Spinfoam Cosmology Setting 

1. Transition Amplitudes from Spinfoams: the General Framework 

Given a boundary graph T and a spin network state \&r on that boundary, a spinfoam model defines possible bulk 
structure as 2-complexes A whose boundary is T and builds a spinfoam probability amplitude [^r] for each of 
these admissible 2-complexes. The boundary graph T and the spin network state define the three-dimensional state 
of geometry and metric on the boundary, while the 2-complcx A and the spinfoam amplitude ,4a defines the bulk 
space-time structure. 

More precisely, the spinfoam amplitude arc defined as local state-sums. One associates algebraic data to the edges 
and faces of A, which usually have an interpretation in terms of discrete space-time geometry. Then all the dynamics 
is assumed to take place at the vertices a of the 2-complex and a local amplitude A„ is defined as a function of the 
algebraic data living on the edges and faces meeting at the vertex a. The spin network state f/'r is understood as 
defining the probability amplitude of the algebraic data on the boundary of A: the data associated to edges (resp. 
faces) of A meeting the boundary will be associated to the vertices (resp. the links) of the boundary graph T. Finally, 
the spinfoam amplitude associated to A is defined as the sum over all possible algebraic data of the product of the 
vertex amplitudes in the bulk and the spin network state on the boundary, which roughly read as: 

-4a t*r] = II ■A<r\jf3<T,ieBo-]ipr(jfedA,ieedA), (65) 

where we have implicitly defined the algebraic data as representations jf of a certain Lie group (usually SU(2) or the 
Lorentz group) on the faces and intertwiners i e between these representations on the edges of A. 

Now considering a fixed graph T, one can use this formalism to define transition amplitudes between two spin 
network states living on that graph. Indeed, introducing the disconnected boundary T U T, one consider a 2-complex 
A whose boundary is T U T and which basically interpolates between an initial copy of T and the final copy of T and 
interprets the spinfoam amplitude as a transition amplitude between a initial spin network state </*{, and a final spin 
network state ipf. 

Going further we will use coherent spin networks peaked on spinor networks as boundary states. These coherent 
spin network states should be considered as coherent wave-packets around the classical spinor networks. Such states 
have been defined in [?} and are based on the coherent intertwiner techniques developed earlier in [l4|, HH, El • Let us 
label such coherent spin network states by their corresponding classical phase space data |{Zg}„ !ee r)- Then spinfoam 
models define transition amplitudes between such coherent spin network states, which can be interpreted as defining 
an effective dynamics for the corresponding spinor networks. This dynamics should correspond to first order at large 
scales to a classical Hamiltonian dynamics plus extra quantum gravity corrections coming from the specific chosen 
spinfoam model for quantum gravity. 
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Our strategy will be to fix the boundary V and the bulk 2-complcx A defining the space-time structure interpolating 
between the initial space slice and the final space slice, and to study the dynamics resulting from the associated 
spinfoam amplitudes. Just the same way, that a fixed simple graph should be enough to describe simple spatial 
geometrical structures and metrics, a simple 2-complex should be enough to describe simple space-time structure 
defining homogeneous dynamics such as the FRW 4-metric. One goal is then to identify which bulk structure A 
corresponds to which mini-superspace model of cosmology. 

A more involved strategy would be to sum over all 2-complexes, as we will comment below from the perspective 
of group field theories for spinfoam models. This would correspond to summing over all possible geometrical and 
topological degrees of freedom for the bulk 4-metric. Here, we are aiming to describing symmetry-reduced physical 
situations which are defined through a finite number of parameters and degrees of freedom. It thus makes sense to 
truncate such a sum over 2-complexes and try to identify the minimal 2-complex faithfully representing the relevant 
physical context and the degrees of freedom of the considered mini-superspace model(s). 



2. Recursion Relations to Differential Equations: to Einstein equations? 

From the procedure described above, we will obtain transition functions ./I^P^im ' z%] depending on the fixed graph 
r, the bulk 2-complex A and the classical initial and final spinor data on the graph T. These transition amplitudes will 
satisfy some differential equations in the spinor variables z's. These differential equations encode both the symmetries 
and the dynamics of the spinfoam amplitudes. Indeed one must keep in mind that symmetries and dynamics are inti- 
mately intertwined in quantum gravity, since the dynamics and evolution are defined by the Hamiltonian constraints 
generating space-time diffeomorphisms. In the context of spinfoam models, this relation is also present. 

More precisely, focusing on the spinfoam model for topological BF theory and restricting ourselves to the gauge 
group SU(2) for the sake of simplicity, one has recursion relations for invariants of SU(2) representations (such as the 
6j-symbol and more generally 3nj-symbols). These recursion relations are understood to be related to the topological 
invariance of the spinfoam amplitudes (see e.g. (l6j). They have also been recently shown to be a quantization of the 
Hamiltonian constraint encoding the evolution and projecting on the Hilbert space of physical states |17| . Furthermore, 
decomposing explicitly the coherent intertwiners and coherent spin network functionals in representations of SU(2), 
one can interpret these coherent spin networks as generating functionals for the SU(2) invariants and corresponding 



spinfoam amplitudes 18[. Then the recursion relations for the the SU(2) invariants and spinfoam amplitudes expressed 
in terms of SU(2) spins get translated into differential equations satisfied by the spinfoam transition amplitudes for 
coherent spin networks on the boundary [T^. These differential equations expresses the invariance of the spinfoam 
amplitudes under certain deformations of the boundary geometry. Therefore they should translate symmetries of the 
considered spinfoam model. But, since some of these deformations induces diffeomorphisms in the time direction, 
they should also encode the dynamics of the theory. 

We will not show this mechanism in details here. We will restrict ourselves to the 2-vertex graph T and to the 
simplest 2-complex A and compute the spinfoam transition amplitudes between homogeneous and isotropic spinor 
data. We will show in this restricted setting that the spinfoam amplitudes does satisfy a differential equation, which 
can be interpreted as a Hamiltonian constraint for homogeneous and isotropic cosmology and from which one can infer 
the classical 4-metric encoded by the spinfoam amplitude. The goal here is to derive the FRW equation and metric 
from this procedure. Applying the same method to generic spinfoam amplitudes, the hope is to derive more general 
differential equations describing the discrete diffcomorphism invariance at the level of 2-complexes and encoding 
general relativity's Einstein equations for 4-metrics. 



3. Spinfoam Transition Amplitudes: a Leading Order Calculation 



In the present paper, we will study states on the 2-vertex graph T with an arbitrary number of edges and choose the 
simplest compatible 2-complex Ai, made of a single vertex in the bulk, following the spinfoam cosmology approach 
introduced in This 2-complcx Ai represents single four-dimensional cell interpolating between the initial 2-vertex 
graph and the final 2-vertex graph, as illustrated in fig. [6] 

As explained above, the spinfoam ansatz in this simplest case of a single spinfoam vertex is that the spinfoam 
amplitude is given by the evaluation of the boundary spin network. Therefore, we choose our spin network states on 
the initial and final 2-vertex graphs and the spinfoam amplitude is the evaluation of these two spin networks. Since 
the boundary is disjoint, we simply get the decoupled product of the evaluations: 

•ASW,^] = #0O^OD- (66) 

In future work, we would like to study the spinfoam amplitudes beyond this leading order 2-complcx with a single 
vertex. Adding vertices and/or using a non-trivial bulk topology, we will then obtain various more complicated 
transition amplitudes beyond this decoupled first order calculation. 
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FIG. 6: Simplest 2-complex Ai interpolating between the initial state Tp* and the final state ^ . Its disconnected boundary 
r U r coincides with the boundary graph of its single vertex a. 



In the next subsection, we will compute explicitly this first order spinfoam amplitude based on Ai for both SU(2) 
BF theory and a spinfoam model for quantum gravity taking into account the simplicity constraints using coherent 
spin networks on the boundary. 



Here, we have taken the point of view of fixing both the boundary graph T on which our spin networks live and the 
bulk spinfoam 2-complex A. Our goal is to compute the corresponding spinfoam amplitudes describing the evolution 
and dynamics of the spin networks for this fixed choice of bulk structure and interpret as a mini-superspace model 
(for cosmology). 

An alternative would be to fix the structure of the boundary but sum over all "admissible" bulks. In order to do 
this, we need to define the list of admissible 2-complexes and to fix their relative weights in the sum. This is done 
automatically by the group field theory formalism which provides us with a non-perturbative definition of the sum 
over spinfoam histories for fixed boundaries (see e.g. [321 133|). 

More precisely, the standard group field theory (GFT) partition function is expanded in Feynman diagrams, which 
are understood as the spinfoam 2-complexes. These are interpreted as the dual 2-skelcton of (pseudo-)triangulations 
made of 4-simplices glued together (along their boundary tetrahedra) , with each 4-simplcx dual to a spinfoam vertex. 
The resulting spinfoam amplitude can be written as a sum over all 2-complexes dual 4d triangulations with fixed 
boundary with a statistical weight: 



where the factor w;[A] is the symmetry factor coming from the Feynman diagram expansion, 7 is the GFT coupling 
constant and #cta counts the number of spinfoam vertices (or 4-simplices) of A. 

Let us revisit our setting with the fixed 2-complex Ai for the 2-vertex graph boundary. It consists in a single 
spinfoam vertex and we have chosen the corresponding spinfoam ansatz for such a configuration given by the straight- 
forward evaluation of the boundary spin network. However, from the GFT point of view, Ai is not dual to a 4d 
triangulation made of a single 4-simplcx. On the contrary, we would need several 4-simplices to get a bulk topologi- 
cally equivalent to Ai. Thus, if we define the spinfoam amplitudes from the group field theory framework, it seems 
that a priori we should not assume the simplest spinfoam ansatz for the amplitude for the 2-complex Ai. On the other 
hand, we should compute the full sum over 4d triangulations compatible with the 2-vertex graph on the boundary 
(more exactly whose boundary is the union of both initial and final 2-vertex graphs). Actually, we do not yet know 
how to control and compute such a sum, which is likely to be divergent, despite the recent progress on this issue. We 
would need to study the coarse-graining of spinfoam amplitudes and the renormalization of the group field theory. 
We would need to extract the relevant interaction terms which should be in the effective group field theory at large 
scale and study the running of the corresponding coupling constants (with the scale of the boundary geometry) . 

One possibility is that the leading order term after renormalization is exactly given by the simplest spinfoam ansatz 
for the 2-complex Ai. This is what happens for the topological BF theory if we assume a trivial bulk topology (which 



4- The Group Field Theory Point of View and the Issue of Renormalization 




(67) 
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seems to dominate the GFT partition function in colored group field theory models (32|) and we can expect this to 
remain true for spinfoam quantum gravity in a low curvature regime. Nevertheless, one should keep in mind the limits 
of our present point of view of fixing the bulk structure and not following the group field theory prescription. 

One ingredient from the GFT perspective which we could keep is the zeroth order term of the expansion. This is 
the identity map between the initial and final boundary coming from the trivial contribution of the 2-complcx directly 
interpolating between the initial and final graph without any vertex. Therefore, we could define a truncated transition 
amplitude as the sum of the identity term plus the first order contribution coming from the one- vertex 2-complex Ai: 

A r A ^[r,^] = A^W^+^^WA'f] = (#|^> +7^(1)^(1)- (68) 

This new zeroth order, with a trivial propagation, gives a totally coupled term between the initial and final boundaries 
while our first order is totally decoupled. In order to focus on the simplest possibility, we will nonetheless not consider 
this zeroth order term but just the first order term when computing the transition amplitude in next sections. 



B. Spinfoam Amplitude and Dynamics for BF Spinfoam 

Let us start applying our program to the spinfoam model for topological BF theory with gauge group SU(2). 
We have already described in great details in section [TT| the classical phase space for SU(2) spin networks on the 
2- vertex graph and its reduction to the homogeneous and isotropic sector through a symmetry reduction by U(iV). 
We have two sets of N spinors, and Wi, attached to the N edges around respectively the vertex a and the vertex 
j3. Both sets of spinors satisfy the closure constraint. Moreover they satisfy the matching constraint along each edge, 
(zi\zi) = (wi\wi). Then the homogeneous and isotropic sector is defined as assuming that the spinors Wi are equal to 
the dual of the spinors z% up to a global phase, \wi] = e 4 * \zi). 

At the quantum level, the components of the spinors are quantized as harmonic oscillators and we recover the 
Hilbcrt space of spin networks on the considered graph (see @, H3, f° r more details). At the end of the 

day, we define coherent spin network states by attaching two coherent intcrtwincrs, and to both the 

source vertex a and the target vertex f3 respectively. These coherent intertwiners are defined as diagonalizing the 
annihilation operators Fij and are labeled by the classical phase space points [ToL [l5j . One can furthermore show 
that they transform covariantly under the U(iV)-action on the spinors and that they provide a decomposition of 
the identify of the space of intertwiners. Finally, they are good semi-classical states, minimally spread about the 
corresponding phase space point. Moreover we understand very well their decomposition in terms of the Livinc- 
Speziale intertwiners (labeled by spins ji € N/2 and spinor variables) introduced in [3l| or in terms of the 

U(iV) coherent intertwiners | J, {z;}) (labeled by the total area J = X^ii an d spinor variables) introduced in (TJ. For 
instance, it will be useful for computational purposes to remind the formula established in (lol [f~5j: 

since the scalar product between the U(iV) coherent intertwiners | J, {z{\) are explicitly known and are homogeneous 
of degree J in the spinor variables [l4| • 

Now, we have our initial 2-vertex graph and our final 2-vertex graph, with the bulk in between defined by the 
2-complcx Ai with a single vertex, as shown in figure [3 Our initial spin network state is the coherent spin network 
labeled with spinors Zi and Wi~. 

i>{z t ,vn}{9i) = ({^ w i}\ ®i 9i\{Zi}), 

where the spin network is a gauge-invariant functional of the SU(2) group elements gi on the edges. Our final spin 
network state is the coherent spin network labeled with spinors and Wi . And the spinfoam amplitude associated to 
Ai is the evaluation of the boundary spin network, which is the decoupled product of the evaluations of the initial 
spin network and of the final spin network: 

AaAzi, Wi,Zi, Wi] = ^. !W .}(I)^ {2 -. iS .}(I) = ({swi}\{zi}) ({swi}\{zi}) , (70) 

where the evaluation of a spin network state on the 2-vertex graph at the identity is simply given by the scalar product 
of the intertwiners living at the two vertices. 
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FIG. 7: Boundary data on the 2-complex Ai defining initial and final boundary coherent spin networks. Focusing e.g. in the 
initial 2- vertex graph, we attach to the vertex a a coherent intertwiner labeled by the spinors {zi} and to the vertex /? another 
coherent intertwiner labeled by the spinors {wi}, both consistently matched to form a coherent spin network ipi z ., w .\. 



Let us focus on a single scalar product, W(zi, iu,) = %ps z . W }(T) = ({swi}\{zi}) . We know how to compute it exactly 

Wfam) = ({Wi}\{zi}) = E J!(./+i)! (J ' { ^ }|J ' {z ' }) 

j 



We can further simplify this expression in the homogeneous and isotropic case by plugging in the relation between 
the two sets of spinors \u)i] = e l ^\zi) and taking into account the closure constraint J^i \ z i)( z i\ = We then obtain 
a formula depending only on the two conjugate variables A and </>: 

w(z -^ ) = ^7!(7Ti)! (Ae " #)2J - (72) 

As mentioned earlier in section Til CI we define the complex variable: 

z = 2V\e~ i %, {A,0} = - => {z,z} =i. (73) 
Then the spinfoam amplitude is simply a Bessel function in z 2 : 

M«,w)=W(z) = EjiTjVl)! (T = ^ h (t) ' (74) 

where I\ is the first modified Bessel function (of the first kind). It is known that Bessel functions satisfy a second 
order differential equations. Here, one can deduce directly from the expansion of W(z) in J that it satisfies: 

C W = 0, with C = z 2 d 2 z + hzd z - z i = d 2 z 2 + zd z - z 4 - 2 = z - 3 d z z 5 d z - z 4 , (75) 

which basically translates the obvious recursion relation on the series' coefficients, (J + 1)!(J + 2)! = ( J + 1)(J + 
2) J!(J + 1)!. 

This constraint C satisfied by the spinfoam amplitude expresses an invariance of the amplitude under certain de- 
formations of the boundary. From a canonical point of view, it should be interpreted as the Hamiltonian constraint. 
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Indeed, it has already been pointed out that the Hamiltonian constraint, as well as the space diffcomorphism con- 
straints, get translated into recursion relations and differential equations satisfied by the amplitudes in the spinfoam 
framework fl6l - [T8T |. We will show below that this is indeed the case and that C is indeed the quantization of the 
Hamiltonian constraint (for BF theory). 

Considering C as the quantization of a classical (Hamiltonian) constraint, let us point out that the two relevant 
terms in the differential operators are z 2 d 2 and z 4 , while the first derivative term zd z and the constant term come 
from ordering ambiguities (for instance, changing the order of the operators in z 2 d 2 clearly affects those terms). 

From this perspective, we also expect this Hamiltonian constraint to express the FRW equation in our present 
cosmological setting. Therefore, we look below at the large scale behavior of the spinfoam amplitude W(z) and of the 
quantum constraint operator C and show their relation to FRW cosmology and to the classical effective dynamics on 
the 2-vertex graph. 



C. Asymptotic Behavior and FRW Equation 

It is well-known that the Bessel function is approximate at large argument by a Gaussian, similarly to a Poisson 
distribution. Let us apply this to our spinfoam amplitude W(z). More precisely, let us perform a stationary point 
approximation of the sum over J. To this purpose, we use the Stirling formula to approximate the factorials for large 
J: 

1 ( Z Y J „ 1 4Jln(f)-2J(lnJ-l) (76) 

J\(J+ 1)1 \2 J J»i 2nJ 2 ' 1 ; 

Then we look for the stationary point(s) of the exponent <p(J) = 4Jln(|) — 2 J (In J — 1) and compute the second 
derivative: 

We get a unique stationary point J = z 2 /4, which justifies that our large J approximation holds when z itself is large 
(in modulus). Then approximating the sum over J by a Gaussian integral around this stationary point, we get: 

w(z) ~ * e 4 

We have checked this formula numerically. What is interesting is that W(z) is simply a Gaussian distribution in z 
with a pre-factor. Under this form, it is obvious that it satisfies a 2nd order differential equation, which is the same 
as W up to a constant shift: 

(C + 3)W(z)=0. (78) 

The fact that it is exactly the same differential operator up to a simple constant shift means that the large scale 
behavior is exactly the same at the classical level and that the deviation from the exact amplitude is at small scales 
in the quantum regime. 

We see that in the asymptotic regime, both the spinfoam amplitude and the constraint operator are very similar to 
the one derived previously in the earlier work on spinfoam cosmology Q, except that we do not use the exact same 

complex variable (z oc \f~\e~ l i here instead of C ~ X — i<f>) and that we easily control the level of approximation that we 
do. Indeed, we know both the exact spinfoam amplitude and constraint operator and their large scale approximations. 
Despite these minor differences, the physical interpretation will be exactly the same. 

Let us focus on the constraint operator C = z 2 d 2 + hzd z — z 4 . Considering the canonical Poisson bracket, {z, z} = i, 
the differential operator d z is the quantization of the classical variable ~z and thus our differential operator turns out 
to be the quantization of the classical Hamiltonian constraint: 

C = z 2 (z 2 -z 2 ) + (5zz). (79) 

Replacing z by its definition 2\f\e~ l i , we see that the first term goes in A 4 and clearly dominates the second term 
which goes in A 2 and can be neglected. This is consistent with the fact that the second term comes from the differential 
operator zd z , which comes from ordering ambiguity in the leading order operator z 2 d z . Therefore, we neglect the 
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term in zz in the definition of our classical Hamiltonian and we factor out the pre-factor z 2 , which leaves us will a 
renormalized classical Hamiltonian derived from the spinfoam amplitude: 

C = z 2 -z 2 = UXsmcf). (80) 

This is exactly the Hamiltonian for FRW cosmology in the simplest case, with vanishing curvature k = and 
cosmological constant A = 0, where we don't have any evolution with <p = and constant A, as we see by comparing 
with the explicit expressions of section Mil More precisely, we have derived the effective dynamics on the 2-vertex 
graph with the Hamiltonian in A sin <j> instead of the usual Xcj> of classical FRW (see section Mil for the comparison 
between the effective and standard Hamiltonians) . 

Going into more details, it is actually to understand how this Hamiltonian constraint comes truly from the dynamics 
of BF theory, as we explain below. Since BF theory is a theory of flat connections, this is the reason that we obtain 
only the flat cosmology case with k = 0. Let us point out that the terms in z~z that we neglected are not with the 
correct scaling in A and do not generate any non-trivial term in k or A. 



D. Recovering the Hamiltonian Constraint 



For now, we have made the link between the spinfoam (transition) amplitude for the quantum dynamics of coherent 
spin networks on the 2-vertex graph and the classical dynamics defined in the earlier section Mil We can actually 
make this link stronger and show that the relation holds at the quantum level and not only in the classical regime 
at large scale A — > oo. More precisely, we show the relation between the differential equation satisfied by the spin 
network evaluation derived in the covariant spinfoam context and the quantization of the classical Hamiltonian on 
the 2-vertex graph defined from a canonical point of view and worked out in Q . 

Considering the evaluation of the coherent spin network state W(zj,iUj) = ip{ Zi/w .y(T), we remind the method 
introduced in [lij and refined in fl7l [l8[ to derive recursion relation or differential equations on the evaluation by 
acting with holonomy operators on the spin network state. Indeed, let us consider the holonomy operator around 
the loop formed by the two edges j and k. It acts by multiplication on the spin network state ip{ ZilWi }(9i)- On the 
other hand, it is compute its action on the labels {zi,Wi} of the spin network state using the recoupling theory of 
representations on SU(2) as discussed in [l6l - [l8j . In short, we get: 

=> 5^W(zi,Wi) = xCffjfl* 1 )^*,™,}^)!^ = 2W(zi, Wi ), (81) 

where the holonomy is taken in the fundamental spin-^ representation (for the sake of simplicity) and 1^9 z r\ is a to- 
be-determined differential operator in the spinor variables. Then by evaluating the action of the holonomy operator 
on the identity, we obtain a differential equation on the spin network evaluation W(zi, w{). 

We can of course consider operators more complicated than a single holonomy operators. However, as soon as it docs 
not act anymore as a multiplication operator in the group elements gi, one has to be careful with the operator ordering 
(e.g. (T^|). but the method still works. Here, we will apply it using a U(7V)-invariant combination of renormalized 
holonomy operator. We consider the classical U(A r )-invariant observable: 

Q = V ( z j\ z j)( z k\zk)(wj\wj)(wk\wk) X(9j9k 1 ) ~ 2 Y V ^ z o\ z o)( z kW){wj\wj){w k \w k ) 
j,k j.k 

= Yl ( 2 ( z j\ z k)(wj\w k ) + [zj\zk)[wj\wk) + (zj\zk](wj\w k ]) - 2 Y^{zj\zj, 

j,k \ j 

= E { 2E iuE} k + F« k Ff k + Fp%) 2 (j2 E%\ (82) 

(83) 

where we used the expression of the group elements gj and g k in terms of the spinor variables in order to compute 
the holonomy X^j^k 1 ) an d where we assumed the matching conditions, (zi\zi)(wi\wi) for all edges i, to simplify the 
constant term. It is clear that this observable vanishes where the group elements are fixed to the identity gi = I: 
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Our strategy will be to compute the action of this operator on the spin network state ip{ Zi , Wi }{9i) at the quantum 
level and deduce the differential equation satisfied by the evaluation W(zj, Wi) by taking g.j = I at the end. 

Let us first see how is the value of Q in the homogeneous and isotropic sector, with \wi] = e 1 ^ for all i's. An 
easy calculation gives: 

Q = 4 A 2 ( 1 + cos 2(j>) -8A 2 = -8A 2 sin 2 ^. (84) 

We recognize both our classical Hamiltonian (|5"4"|) on the 2-vertex graph for the special values 7° = 7 and 7 1 = 
and the square of the constraint C derived in (|80[) from the spinfoam amplitude. This makes the link between the 
canonical formalism with a Hamiltonian and the covariant perspective with the spinfoam amplitude satisfying certain 
constraint. 

Let us now work out the quantized version Q and check that its action on fc jW s(ft) vanishes on trivial holonomics 
gi = I. This will also provide us with the exact differential equation gcnerically satisfied by the evaluation W(zj, lUj), 
even without assuming the homogeneous and isotropic ansatz. We follow the quantization procedure for the spinor 
variables: holomorphic coordinates Zi and Wi are quantized as multiplication operators while the anti-holomorphic 
variables z 2 ; and Wi respectively become the differential operators d Zi and d Wi . This leads to the following quantization 
for the E and F observables (7l.Hoj: 

E% = zld z o+zld z]l 

Fjk = ( z< j z k ~ z j z k)> 

F 3 V = (d z0 d zl -d z] d z o), (85) 

and similarly for the operators acting at the vertex f3 as differential operators in the lUj's. We can then compute the 
action of these operators on W(zi,Wi) = 7 (dct X) J / J!( J + 1)! where X is the following 2x2 matrix: 



i jk 

Following the natural notation as in 0, we will denote the functional (detX) J as the quantum state \J). Then the 
three components of the Q operator, E" k Ej k , F" k Fj k and F^F^ , respectively act as a number of quanta operator, 
a creation operator and an annihilation operator in the \J) basis. More precisely, after a lengthy but straightforward 
calculation, we get: 

Y,E? k El k \J) = 2J(J + N-2)\J) (86) 
j,k 

T. F ^ F "k\J) = 2|J + 1) 

j,k 

J2 P ^ P Jk\J) = 2J(J + l)(N + J-l)(N + J-2)\J-l). 

It is worth pointing out that the same expressions were computed earlier in a more elegant way in Q by working out 
the commutation relations between these operators in the U(iV)-invariant space. Putting these results together, we 
define the quantized version of the Q observable: 

Q ee £ (2E? k Ef k + tf k Ff k + FtfFfc) - 2{E + N — l) 2 - 2(N - 1), (87) 

where E = J2i E?i is shown to act as E\J) = 2J\J). The last term (E + N - l) 2 + (N - 1) corresponds to the 
quantization of the classical term . E^) 2 up to sub-leading contributions interpreted as ordering ambiguities. This 
specific choice gives the expected result: 

ew(^^) = QEj0fij! =«■ ( 88 ) 

This is a fourth order differential equation on the spinfoam amplitude yV(zi, Wi), which correspond to the quantization 
of the classical Hamiltonian on the 2-vertex graph. 

It is further possible to write all the differential operators in terms of derivative of detX, and thus in terms of the 
single complex variable z, then the operator Q should give the square of the differential equation C 2 (up to ordering 
terms) satisfied by the spinfoam amplitude. 
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E. How to Depart from Flat Cosmology? 



Working with the spinfoam model for topological BF theory with gauge group SU(2), we have derived the differential 
equation satisfied by the spinfoam amplitude for fixed boundary (the 2-vertex graph T) and fixed bulk (the 2-complex 
Ai) and shown its relation to the classical Hamiltonian for the dynamics on the 2-vertcx graph. The Hamiltonian that 
we have obtained corresponds to the flat FRW cosmology (in vacuum) . The natural question is whether it is possible 
or not to obtain the models of effective dynamics of section IIHI for non- vanishing curvature k ^ and cosmological 
constant A ^ with this spinfoam cosmology framework. 

From the previous analysis, we understand that the flatness of the Hamiltonian k = = A, or equivalently 7 = 7 
and 7 1 = in the cosmological Hamiltonian ansatz (|54[) . comes from the definition of the spinfoam amplitude as the 
evaluation of the boundary spin network at the identity. This comes from working with the spinfoam path integral 
for topological BF theory, which projects onto trivial connections and holonomies. Therefore, in order to obtain 
non-flat FRW cosmology, it is clear that we have to depart from BF theory and the trivial spinfoam ansatz defining 
the amplitude as the evaluation of the boundary spin network state at the identity. For instance, it is natural to 
expect that evaluating the boundary spin network on non-trivial holonomies will produce curvature and lead to FRW 
cosmology models with non- vanishing curvature and cosmological constant. 

In order to derive such a new spinfoam amplitude, we see two non-exclusive possibilities. First, we can change the 
spinfoam model. Either, we can attempt to modify the spinfoam amplitudes by hand and introduce curvature and 
cosmological constant terms in the BF amplitude. Or we can use a spinfoam model built for quantum gravity, such as 
the EPRL-FK spinfoam amplitudes [3 US H3] ■ However, if we use such a model, the standard spinfoam amplitude 
ansatz is still to evaluate the boundary spin network around each vertex at the identity. Therefore, we expect no 
difference from topological BF theory if our bulk 2-complex contains a single vertex as considered here with our 
simplest 2-complex Ai. We illustrate this by working out in the next subsection the spinfoam transition amplitude 
for the 2-vertex graph still with the single- vertex 2-complex Ai for the Dupuis-Livine variant of the EPRL model 
based on the holomorphic simplicity constraints [1(1 [H| • Then the second alternative is to change the bulk 2-complex, 
and even the boundary graph, in order to allow for more intricate space-time structure with non-trivial topology and 
geometry By allowing more than one vertex in the bulk 2-complex, we expect to have curvature generated through 
the gluing of the dual 4-cells. 

Here, we will show how to modify the spinfoam amplitude and depart from the mere evaluation of the boundary spin 
network state in order to obtain non-flat FRW cosmology. We will use the realization of the Hamiltonian constraint as 
a differential equation satisfied by the spinfoam amplitude. Indeed, we have derived the differential equation satisfied 
by the evaluation of the boundary spin network state at the identity. We can now modify this differential equation 
to take into account non-vanishing curvature and cosmological constant and investigate how the spinfoam amplitude 
changes. 

Calling e ee 2 £\ k Ef k E^ k - 2(E + N - l) 2 - 2(N - 1), / = Ff k F^ k and ft = F^F^ the operators introduced 
above and defined as differential operators in the spinor variables Zj and Wi, the differential equation satisfied by the 
coherent spin network evaluation was simply defined by the operator Q = e + / + ft: 

Q ( W|{zJ) = Q £ J[{ j +iy (dot X) J = . (89) 

Following our analysis of the classical dynamics on the 2-vertex graph in section IIII1 we consider the modified differ- 
ential operator = e + n(f + ft) for k G R. The corresponding classical constraint (in the large A regime) in the 
homogeneous and isotropic sector is: 

q(k) = _ 2 2A 2^ _ K cos ( 90 ) 

which corresponds to our effective FRW Hamiltonian (|55|) with parameters 7 = 1, 7 = k and 7 1 = (up to a global 
(— 2)-factor), thus leading to modified FRW cosmology with non-vanishing curvature k 7^ 0. We can easily take into 
account a non- vanishing cosmological constant by adding a E 3 term in our constraint. But for the sake of simplicity, 
we will simply describe here the case 7 1 = 0. 

Let us then solve the constraint QS K > — at the quantum level. Noting as before |J) ee (dctAT)' 7 , we apply the 
modified operator Q,( K > to a linear combination of the basis vectors | J) and we obtain a 2nd order recursion relation 
on the coefficients: 



a (K) 

I J) = <=> 



y J\(J + i)\ 



(J + N)(J + N- l)a { j% - I [2 J( J + N) + N(N - l)}a^ } + J{J + l)a^x = 
«(JV + l)ai = (iV-l)ao 

(91) 
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The spectral properties of similar operators were studied in [6| when analyzing the Hamiltonian dynamics of the 
quantum 2- vertex model. Applying the same techniques, we look at the asymptotics of the recursion relation at large 
J and solve it at second order (in 1/J). We then get: 

a J ^ T ,l ; cos a; J, coscj=i. (92) 

Wc have checked the accuracy of this asymptotics numerically using Maple 15. For k = 1, the oscillation frequency u 
vanishes and we recover the flat case which we have already described. For \n\ < 1. the frequency uj becomes purely 
imaginary and we have an exponential solution instead of the oscillating behavior. 

We can modify slightly the action of the operators e, /, / ' on the states | J) in order to obtain an exact analytical 
expression for the physical state. For instance, we can take: 

e|J) = -4J 2 |J) 

/| J) = 2\J + 1) (93) 
/t|J) = 2J 2 (J-1) 2 |J-1) 

o(«0 - -. 

and solve the modified equation Q = e + n{f + / T ) =0. Then we get the states: 

v-^ coswJ . . ^( K ) v-^ coswJ _ , , , 

2 E^7i^l J ) = 2 E^7i^( dctx ) =°- ( 94 ) 

-:(«) 

This new operator Q has the same classical limit than our original operator but can be interpreted as 

differing in operator ordering. It is simple to define these modified operators e, /, /' as differential operators in the 
single complex variable detX, or equivalcntly z, but we lose the direct translation as differential operators in the 
original spinor variables which allow to define these operators as acting on the whole phase space and not only in the 

~ -+ ~ M 

homogeneous and isotropic sector. Moreover, this hides the relation of the operators e, /, / ' , Q with the holonomy 
operators and the interpretation of the solution state as a coherent spin network evaluation. 

Even if we do not have an explicit closed formula for the coefficients a j , this gives the spinfoam amplitude than 
we should define in order to get non-flat FRW cosmology. The next step is to understand how this amplitude can be 
defined as the evaluation of the coherent spin network on non-trivial holonomies or with an operator insertion, i.e of 
the type ({^WiWO^ \{zi}}. Already the asymptotics of the coefficients ctj provide us with some clues about how to 
realize this. Indeed, the series based on the asymptotics can be easily realized as the evaluation of the coherent spin 
network on some non-trivial holonomies: 

= \ ^ e iu \{ Zi }) + ( W| ®, e~^\{ Zl })] (95) 

where we act with the same U(l) group element e lul , resp. e~ lu , on all the legs of the 2-vertex graph. This corresponds 
to the evaluation of the coherent spin network states on SU(2) group elements defined as the 2x2 diagonal matrix 
[e lu , e~ lu ] (resp. [e~ lu ,e +lu ]) in the (|z.;), basis on each edge of the graph. As we have shown above, this 
functional satisfies a differential equation, whose classical counterpart at large scale, is and thus corresponds to 

our Hamiltonian (|55p for the modified FRW cosmology with non-vanishing curvature. 

It seems possible to follow the same strategy to take into account a cosmological constant. We expect to find a 
spinfoam amplitude similar the one postulated in Q or/and related to the g-deformation of the spinfoam amplitudes. 
Then the goal will be to interpret it as some particular evaluation or operator insertion of the boundary coherent spin 
network state. We postpone such a detailed analysis to future investigation. 



F. Cosmological Dynamics with Holomorphic Simplicity Constraints 

Let us apply our method to analyze the spinfoam amplitude for the Dupuis-Livinc (DL) spinfoam model for 
Ricmannian quantum gravity based on holomorphic simplicity constraints. We will not review the details and the 
definition of this model, but we will refer the interested reader to [l(| for the full presentation or to [38| for a quick 
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summary of the model's basic features. We will only point out here that the resulting spinfoam amplitudes are very 
similar to the ones of the EPRL-FK spinfoam models j35l . |37| , except that the diagonal simplicity constraints are not 
strongly imposed and that the wave-packets for the boundary states are thus slightly enlarged Gaussian distributions 
compared to the coherent states used for the EPRL-FK models. This allows for an easier analysis and understanding 
of the symmetries, constraints and amplitudes of the model. 

The DL model is based on the gauge group Spin(4) = SU(2)l x SU(2)a; made of two copies of SU(2). The classical 
phase space on a fixed graph is therefore simply two copies of the SU(2) spinor network phase space and the boundary 
spin networks at the quantum level are simply tensor products of spin network states for the left copy of SU(2) and the 
right copy. The the holomorphic simplicity constraints are easily imposed at both the classical level and the quantum 
level. The resulting structure is coherent spin network states labeled by two sets of spinors living on the graph, the 
z ' R, s respectively from the left and right sectors: 

*{»*-,«* S }(Se»fl?) = 1>{z*-Z}(9e)1>{,*i}(g?), (96) 

where {g L ,g R ) € SU(2)^ x SU(2)# represents a Spin(4) group element. Then the holomorphic simplicity constraints 
simply impose that all the right spinors are equal to their left counterpart up to a fixed factor, Vv,e, z Rv e = pz Lv e 
where p > is related to the Immirzi parameter [Ioj |. Therefore, removing the useless subscript L or R for the spinor 
variables, the DL spinfoam model is based on the coherent spin network states for Spin(4) defined as: 

*{*« e }(9e,9?) = ^}G£)Vw(0f). (97) 

Repeating our derivation of the spinfoam amplitude for the 2- vertex graph and Ai 2-complex in the bulk, we obtain 
once more a decoupled spinfoam amplitude which is the product of the initial and final spin network evaluations: 

A^[zi,Wi,%,Wi\ = W DL (zi,Wi)W DL (fr t Wi) t 



with W Uh {z uWi ) = * {Zi , Wi} (l) = ip {ZitWi} (I)ip {pZi , pWi} (I) = ({WiWiziMtepWiMpzi}) ■ (98) 

Since these scalar products are known analytically, we can compute exactly these amplitudes. Nevertheless, it is 
simpler to have a look at the asymptotic behavior. Using the stationary point approximation derived earlier (|77[) , we 
get for the homogeneous and isotropic sector: 

W^(z) ~*e4 e ^ = * e U+P 2 )4 . (99) 

Z'^y/lT {pZj^y/TT Z p a TT 

This Gaussian wave-packet is once more solution to a 2nd order differential equation whose leading order is of the 
type d\ — z 2 (up to some global power of z). Considering the large scale regime, we obtain a classical constraint 
identical to the BF theory case corresponding to the flat FRW cosmology with no dependence (at leading order) on 
the Immirzi parameter p. 

This confirms that we should study more complicated bulk configurations or change the spinfoam amplitude ansatz 
for the simple single- vertex bulk Ai in order to obtain some effective dynamics for more generic cosmology, with 
possibly curvature and cosmological constant. 



V. CONCLUSION AND OUTLOOK 



We have discussed the general procedure of truncating loop quantum gravity and spinfoams to a fixed finite graph 
and studying the dynamics of the geometry on that fixed graph at both classical and quantum levels. The goal is to 
define mini-superspace models that can be interpreted as cosmological models or as toy models to investigate and test 
the spinfoam dynamics. Following this logic, we have studied in detail the kinematics and dynamics of loop gravity 
on a 2-vcrtex graph. We have defined the dynamics from both the canonical point of view and using the covariant 
spinfoam amplitudes and showed the consistence between the two approaches. We have further defined the reduction 
of a homogeneous and isotropic sector at the kinematical level and showed that the dynamics of this sector can be 
understood as some modified effective Friedmann-Robertson- Walker (FRW) cosmology. This confirms and improves 
the earlier results about the 2- vertex model derived in 0, 0-0] . These results show that the present approach might be 
very promising to derive mini-superspace models from spinfoam models and to study the quantum gravity dynamics 
and fluctuations in these restricted settings. 

Concerning the interpretation of the 2- vertex model as a cosmological mini-superspace model, the limitations of our 
approach are clear. First, we absolutely need to include matter in our analysis, both at the canonical level (following 
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the earlier work Q) and in the spinfoam calculations, in order to obtain a realistic FRW model. Second, we have 
identified and explained some discrepancies in the large scale behavior of our 2-vertex dynamics with the standard 
FRW dynamics when we turn on the curvature k or the cosmological constant A. This shows the limits of our current 
approach. 

Moreover, concerning the derivation of the 2-vertex cosmology from spinfoam amplitudes, we also see two main 
issues. On the one hand, we have discussed the effects that curvature and cosmological constant cause on the spinfoam 
transition amplitude, but we still need to understand how to encode them in the coherent spin networks living in the 
boundary of the two-complex. On the other hand, we need to study the transition amplitudes for more complicated 
bulk structures, beyond our calculations for a simple one 4-cell bulk. This involves studying more the spinfoam 
amplitudes and their coarse-graining/renormalization in order to understand how the transition amplitudes for more 
complex bulks look like and to determine if the one 4-cell bulk is actually the dominating contribution. Both issues are 
intertwined since we expect that considering more complicated bulks would (at least) induce curvature. Actually, for 
the simplest 2-complex that we have considered, the boundary of the spinfoam vertex coincides with the boundary of 
the 2-complex, and then the transition amplitude is just the evaluation of the boundary spin network on the identity 
(according with the spinfoam ansazt), meaning that we can only describe a flat model (BF theory) in this simple 
setting. Nevertheless, besides these questions to solve, we have shown how to derive differential equations satisfied by 
spinfoam amplitudes and how they are equivalent to the Hamiltonian constraint, in the context of the 2- vertex graph. 
We hope that this method will be generalizable to more complicated graphs and bulk structures. 

Beyond these issues to solve, we would like to point out three promising directions of work. First, we would like 
to investigate the freedom in modifying the one 4-cell spinfoam amplitude and study the resulting dynamics and 
corresponding Hamiltonian. Second, we would like to apply our methods to more complicated graphs, for example 
the 3+N-vertex graph proposed in Q or possibly graphs with an infinite number of vertices, in order to test the 
validity of truncating loop (quantum) gravity to a fixed graph beyond the simplest 2-vertex graphs and to try to 
generate more realistic cosmological models, beyond the homogeneous and isotropic FRW cosmology. For instance, 
we hope to identify some simple family of graphs with a number of vertices that can be sent to infinity in order to 
model midi-superspace models for cosmology with inhomogeneities. Finally, it could be interesting to develop some 
simplified group field theory for our 2-vcrtcx graph model, which would admit specific boundary graph and would 
sum over a restricted set of 2-complexes, and compare it to the group field theory for the EPRL-FK spinfoam model 
[39j and to the spinfoam amplitudes for loop quantum cosmology poMi^ . 
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Appendix A: Spinors and Notations 

In this appendix, we introduce spinors and the related useful notations, following the previous works [f| [l5j . 
Considering a spinor z, 



(z\ = ( z° z 1 ) , 



we associate to it a geometrical 3-vector V(z), defined from the projection of the 2x2 matrix \z)(z\ onto Pauli 
matrices a a (taken Hermitian and normalized so that (er a ) 2 = I): 

\z){z\ = U(z\z)I + V(z)-a). (Al) 



2 

The norm of this vector is obviously |V(z)| = (z\z) = \z°\ 2 + \z 1 \ 2 and its components are given explicitly as: 

y* = |z°| 2 -|*T, V x =2 5R(z°z 1 ), V y =2S(z°z 1 ). (A2) 
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The spinor z is entirely determined by the corresponding 3- vector V(z) up to a global phase. We can give the reverse 
map: 



where e 1 ^ is an arbitrary phase. 

Following 0, we also introduce the map duality <; acting on spinors: 



M i = Jo > ^ = - L ( A4 ) 



a; 



This is an anti-unitary map, (^z\<^w) = (w\z) = (z\w), and we will write the related state as 



\z\ = s\z), [z\w\ = {z\w). 
This map c maps the 3-vector V(z) onto its opposite: 

\z][z\ = \((z\z)I-V{z)-a). (A5) 

Finally considering the setting necessary to describe intertwiners with N legs, we consider N spinors z% and their 
corresponding 3- vectors V(zi). Typically, we can require that the N spinors satisfy a closure condition, i.e that the 
sum of the corresponding 3-vectors vanishes, ^ V( z i) = 0. Coming back to the definition of the 3- vectors V(zi), the 
closure condition is easily translated in terms of 2 x 2 matrices: 

J2\z i ){z i \ = A(z)l, with A(z) = ^{z i \z i ) = ^\9{z i )\. (A6) 

i i i 

This further translates into quadratic constraints on the spinors: 

E^°4 = 0, ^|^| 2 = ^K| 2 = ^). (A7) 

i i i 

In simple terms, it means that the two components of the spinors, zf and z\, are orthogonal N- vectors of equal norm. 



Appendix B: Effective Hamiltonian of the FRW model in LQC 

In LQC, following LQG, one starts by describing the metric in triadric form. Therefore we introduce a fiducial 
triad being i an internal SU(2) index, and write the fiducial metric of uniform curvature as °q a b = "w^w^y. 
Actually, °u> l a denote a basis of left-invariant one forms that satisfy the Maurer-Cartan equation 

d a °ul - —C^M^l (Bl) 

The three different spaces, defined via Eq. P5|) . can alternatively be described with the following conventions: 

• Euclidean case: °ui l a = 8 l a and = Sij. 

• Spherical case: In this case kij is the Cartan-killing metric of su(2) and °q a b is the metric of the round 3-sphere 
of radius 2. 

• Hyperbolic case: It is characterized by having structure constants Cj fc = S^Ski — 8 l k Sji, and kij = 5^. 

As in LQG, in LQC the phase space is described in terms of the densitized triad, Ef, and the Ashtekar-Barbero 
connection, A % a = T l a +jK^, being the spin connection, 7 the Immirzi parameter, and K % a the extrinsic curvature 
in triadric form. For the FRW model under consideration the densitized triad is parameterized in terms of a single 
parameter p, in the form Ef = pV 2 ^ 3 \/°q°ef. Here °ef denotes the fiducial co-triad, such that = S^Sf. On 

the other hand, we parameterize the Ashtekar-Barbero connection in the following way 
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cV ~ 1/3 °GJ l a . 



• Euclidean case: A l a 

• Spherical case: A l a = + J Vo 

• Hyperbolic case: A l a = — ej l °luI + cV 1 ^ 3 °oj t a , (e£: totally antisymmetric symbol). 

— 1/3 

This parameterization is such that, in all the cases, = cV "w^, being °uj l a the fiducial triad corresponding 
to each case. Also, according with the above definitions, we have {c,p} = 8irG/3. For more details we refer the 
interested reader to e.g. ^JtfXl ■ 

The canonical transformation from the variables (a, n a ) to the variables (p, c) is given by 



4^G 7 y o 2/ \/' 



(B2) 



where r Q = 1 in the flat and open models, and r Q = 1/2 in the closed model. 

The effective Hamiltonian derived in LQC, in terms of some dimensionless quantity \i whose meaning we will discuss 
later, is given by 



• Flat case \M \M 4 



k = 0, r a = 1 



N 



16ttG 



6 ^sin 2 (uc) „, o/o 
— VP f^ + 2Ap 3 / 2 



(B3) 



• Open case [44j, |45|: k = -1, r a = 1 



• Closed case W6l, m 



1, r, 



N 



H, 



eff 



2 3 / 2 A 
16ttG 



16ttG 

= 1/2 
G 



2 2 



v ^[sin 2 (H- 7 VK 2/3 



7 2 /i 2 



VP 



sin 2 (/ic) — sin 2 



1/3 ^ 



2Ap 



3/2 



+ (l+7 2 ) 



A 



3/2 



(B4) 



(B5) 



Note that in the limit /i -> wc have H c g — > Hfrw with 



H 



N 



FRW 



6 



[c 2 + k^V^rl 



2Ar> 



3„3/2 



where we have done the approximations 



sin 



'(lie) 



sin 2 (/^/ 1/3 /2) Vj 



2/3 



(B6) 



(B7) 



Introducing the canonical transformation (|B2[) in Eq. (|50|) , it is straightforward to check that the above Hamiltonian 
is indeed the Hamiltonain of the FRW model in the variables (c,p). Therefore the above effective Hamiltonians 
introduce corrections to the standard dynamics of the FRW model. Those corrections come from the quantization 
and are unimportant as long as the parameter /i is small. 

Actually, in the literature of LQC different values of /j, have been considered. In the first analyses of the dynamics 
of the FRW model /i was fixed to have a constant value fi a [27], [28| . Later on the dynamics was "improved" by doing 
this parameter dynamical and given by jl oc p" 1 / 2 [29]. In this way, the limit Ji — >• corresponds to the limit of large 
areas, and then the effective "improved" dynamics coincides with the standard one when the cell of the universe under 
consideration is sufficiently large. 

Let us now translate the effective Hamiltonians of LQC into our variables (A, </>). From the canonical transformations 
given in Eqs. ([56]) and (|B2[) we readily obtain the relations between the variables (c,p) and the variables (A, </>): 



p = — A, 



372 ' 



(B8) 



Introducing this transformation and choosing again N = A 3 / 2 we get 
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• k = 0: 



I 



6 2 sin 2 (/icr</>) 



r 



2AZ 2 A 3 



(B9) 



fc = -1: 



/ 

16ttG 



sin 2 (^)- 7 VV 2 / 3 



2AZ 2 A 3 



(BIO) 



jfc = 1: 



47rG 



7 2 /j, 2 



sin 2 (n<J<p) — sin 2 



1/3 1 



(i+r) 



2/3 




(Bll) 



Now, in the case of the "improved dynamics", the parameter \x is given by Jx oc A -1 / 2 . 

Note that fx needs to be constant for the above Hamiltonians to admit a natural implementation in the spinor 
framework, in the sense that they can be writen as a function of the area A and of the basic holonomy elements of 
the spinor framework, h$ = e 2 ^, In fact, in that case we can write 

hf - hT^" 

sm(^) = * . (B12) 
li 

On the other hand, the effective Hamiltonian of the "improved dynamics" does not admit a natural implementation in 
our framework. Indeed the expression sin(/kr<^>) with /2 oc A -1 / 2 is not in the algebra generated by the basic variables 
A and !%&. Nonetheless, the combination /2cr</> turns out to be canonically conjugate to the volume (up to a constant 
factor) [29[. Therefore, the "improved" Hamiltonian would be naturally implementable within the spinor framework 
if, after the symmetry reduction, the overall variable were the volume instead of the area. While the 2- vertex graph 
that we are considering in this paper is not suitable to hold the improved dynamics of LQC, it might be that a more 
complicated graph and a suitable implementation of the symmetry reduction would indeed pick up the volume as the 
overall variable. We leave the investigation of this possibility for the future. 
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